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Abstract 

This report provides the results of an extensive and important study of the potential for a U.S. scientific 
program that will extend our knowledge of neutrino oscillations well beyond what can be anticipated from 
ongoing and planned experiments worldwide. The program examined here has the potential to provide the 
U.S. particle physics community with world leading experimental capability in this intensely interesting and 
active field of fundamental research. Furthermore, this capability is not likely to be challenged anywhere 
else in the world for at least two decades into the future. The present study was initially commissioned 
in April 2006 by top research officers of Brookhaven National Laboratory and Fermi National Accelerator 
Laboratory and, as the study evolved, it also provides responses to questions formulated and addressed to 
the study group by the Neutrino Scientific Advisory Committee (NuSAG) of the U.S. DOE and NSF The 
participants in the study, its Charge and history, plus the study results and conclusions are provided in this 
report and its appendices. A summary of the conclusions is provided in the Executive Summary. 
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1. EXECUTIVE SUMMARY 



This report provides the resuhs of an extensive and important study of the potential for a U.S. 
scientific program that will extend our knowledge of neutrino oscillations well beyond what can 
be anticipated from ongoing and planned experiments worldwide. The program examined here 
has the potential to provide the U.S. particle physics community with world leading experimental 
capability in this intensely interesting and active field of fundamental research. Furthermore, this 
capability is not likely to be challenged anywhere else in the world for at least two decades into 
the future. The present study was initially commissioned in April 2006 by top research officers of 
Brookhaven National Laboratory and Fermilab and, as the study evolved, it also provides responses 
to questions formulated and addressed to the study group by the Neutrino Scientific Advisory 
Committee (NuSAG) of the U.S. DOE and NSF. The participants in the study, its Charge and 
history, plus the study results and conclusions are provided in this report and its appendices. A 
summary of the conclusions is provided in this Executive Summary. 

The study of neutrino oscillations has grown continuously as its key impact on particle physics 
and various aspects of cosmology have become increasingly clear. The importance of this fun- 
damental physics was recognized by the National Research Counciiyj and the Office of Science 
and Technology Policy [|21, and its national budget priority has been established in a joint OSTP- 
0MB policy memorandum in 2005 (Si. In fact, as the present study confirms, it is now possible 
to design practical experiments that are capable of measuring all the parameters that characterize 
3-generation neutrino oscillations, including the demonstration of CP-violation for a significant 
range of parameter values beyond present limits. Also, one of the experimental approaches, in 
which the detector (regardless of technology) is deployed deep underground, considered in this 
study has the potential to contribute, to a significant improvement of our knowledge about nucleon 
decay and natural sources of neutrinos. 

The two experimental approaches studied here are complex in their detailed technical realiza- 
tion, comprising several detector technologies, various specific neutrino beam designs and different 
measurement strategies. They have in common, however, the exploitation of experimental base- 
Unes of ~1000 km (a key advantage of a U.S. based program) and both approaches make effective 
use of existing Fermilab accelerator infrastructure with modest upgrades. The experimental de- 
tectors required are very massive (in the several hundred kiloton range) because the interaction 
rates are small. The designs for such detectors vary from already-demonstrated at a scale of 50 
kTon (Super Kamiokande) to somewhat speculative (large liquid Argon). In both cases, signifi- 
cant R&D is still needed to demonstrate feasibility and obtain a reliable cost estimate for the scale 
needed here. The study has shown however, that it will be feasible and practical to carry out the 
desired program of important neutrino physics, perhaps together with improved nucleon decay and 
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natural neutrino investigations in the same neutrino detector. 

The output of the present study is twofold: 1) technical results and conclusions that report the 
results of the study and address the charge letter; 2) answers to the 15 questions posed to the study 
group by the NuSAG Committee. These two outputs comprise more than 50 pages of detailed 
commentary and they are provided in full in the body of the report and Appendix A. Here, we 
attempt to provide a somewhat condensed version of the study results and conclusions while urging 
the reader to consult the full text of the report on any points that may appear to be questionable or 
unclear. The summary results and conclusions were discussed and agreed to at the September 17, 
2006 meeting of the study group. 

Results and Conclusions: 

• Very massive detectors with efficient fiducial mass of > 100 kTon are needed for the ac- 
celerator long baseline neutrino program of the future. We define efficient fiducial mass as 
fiducial mass multiplied by the signal efficiency. For accelerator based neutrino physics, this 
could correspond to several hundred kTon if the detector is a water Cherenkov detector and 
> 100 kTon if it is liquid argon TPC with high expected efficiency. These detectors could 
be key shared research facilities for the future particle, nuclear and astrophysics research 
programs. Such a detector(s) could be used with a long baseline neutrino beam from an 
accelerator laboratory to determine (or bound) leptonic CP violation and measure all param- 
eters of neutrino oscillations. At the same time, if located in a low background underground 
environment, it would have additional physics capabilities for proton decay and continuous 
observation of natural sources of neutrinos such as supernova or other astrophysical sources 
of neutrinos. 

• The Phase-II program will need considerable upgrade to the current accelerator intensity 
from FNAL. Main Injector accelerator intensity upgrade to ~ 700 kW is already planned 
for Phase-I of the program (NOvA). A further upgrade to 1.2 MW is under design and 
discussion as described briefly in this report. The phase-II program could be carried out 
with these planned upgrades. Any further improvements, perhaps with a new intense source 
of protons, will obviously increase the statistical sensitivity and measurement precision. 

• A water Cherenkov detector of multi-lOOkTon size is needed to obtain sufficient statistical 
power to reach good sensitivity to CP violation. This requirement is independent of whether 
one uses the off-axis technique or the broadband technique in which the detector is housed 
in one of the DUSEL sites. 

• High signal efficiency at high energies and excellent background reduction in a liquid argon 
TPC allows the size of such a detector to be smaller by a factor of 3 compared to a water 
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Cherenkov detector for equal sensitivity. Such a detector is still quite large. 

• The water Cherenkov technology is well established. The issues of signal extraction and 
background reduction were discussed and documented at length in this study. The needed 
background reduction is achievable and well understood for the broadband beam discussed 
in this report, but not yet fully optimized. Key issues for scaling up the current generation 
of water Cherenkov detectors (Super- Kamiokande, SNO, etc.) and locating such detectors 
in underground locations in DUSEL are well understood. The cost and schedule for such a 
detector could be created with high degree of confidence. A first approximation for this was 
reported to the workshop. 

• For a very large liquid argon time projection detector key technical issues have been iden- 
tified for the building of the detector. A possible development path includes understanding 
argon purity in large industrial tanks, mechanical and electronics issues associated with long 
wires, and construction of at least one prototype in the mass range of 1 kTon. 

• In the course of this study, we have examined the surface operation of the proposed mas- 
sive detectors for accelerator neutrino physics. Water Cherenkov detectors are suitable for 
deep underground locations only. Surface or near-surface operation of liquid argon TPCs 
is possible but requires that adequate rejection of cosmic rays be demonstrated. Surface or 
near-surface operation capability is essential for the off-axis program based on the existing 
NuMI beam-line because of the geographic area through which the beam travels. 

• Additional detailed technical conclusions of the study are noted in the Results and Conclu- 
sions section of this report. These results could influence the detailed design of the specific 
program selected. 

Detailed sensitivity estimates for the choices under consideration can be obtained from Section 



13 Here we will give a broad comparison of the different experimental approaches. 

In the course of this year long study we have been able to draw several very clear conclusions. 
Regardless of which options evolve into a future program, the following will be required. 

1 . A proton source capable of delivering 1-2 MW to the neutrino production target. 

2. Neutrino beam devices (targets and focusing horns) capable of efficient operation at high 
intensity. 

3. Neutrino beam enclosures which provide the required level of environmental and personnel 
radiological protection. 
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4. Massive (>>100 kton) detectors which have have high efficiency, resolution and back- 
ground rejection. 

5. For each of the above items, significant investment in R and/or D is required and needs to be 
an important aspect of the current program. 

We have found that the main areas of this study can be discussed relatively simply if we divide 
them into two broad categories : 1) The neutrino beam configuration and 2) The detector technol- 
ogy. Further, we are able to summarize our conclusions in two tables which show the pros and 
cons of the various options. 

In Table|I]we compare the pros and cons of using the existing NuMI beam and locating detectors 
at various locations, versus a new wide band neutrino beam, from Fermilab but directed to a new 
laboratory located at one of the potential DUSEL sites, i.e. at a baseline of 1300 to 2600 km. 

In Table|ll]we compare the pros and cons of constructing massive detectors ( 100 - 300 kT total 
fiducial mass) using either water Cherenkov or liquid argon technology. 
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Pro 


Con 


NuMI On-axis 


Beam exists; 
Tunable spectrum; 


L ~ 735 km 

Sensitivity to mass hierarchy is limited 
Difficult to get flux < 3 GeV 


NuMI Off-axis 
(1st maximum) 


Beam exists ; 
Optimized energy; 
Optimized location for 
1st detector; 

Site will exist from NOvA project; 


L ~ 800 km 

Limited sensitivity to mass hierarchy 


NuMI Off-axis 
(2nd maximum) 


Beam exists; 
Optimized energy; 
Improves mass hierarchy 
sensitivity if 613 is large; 


L ~ 700-800 km; 
Extremely low event rate; 
A new site is needed; 
Energy of events is ^ 500MeV; 
Spectrum is very narrow 


WBB to DUSEL 


More optimum (longer) baseline; 
Can fit oscillation parameters 
using energy spectrum; 
Underground DUSEL site for detector; 
Detector can be multi-purpose; 


New beam construction project >$100M; 
Multi-year beam construction; 



TABLE I: Comparison of the existing NuMI beam to a possible new wide band low energy (WBLE) beam 
to DUSEL 
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Pro 


Con 


Water 


Well understood and proven technology; 


Must operate underground; 


Cherenkov 


Technique demonstrated by SuperK (50kT); 


Scale up factor is < 10; 






Cavern stability must be assured 






and could add cost uncertainty; 




New background rejection techniques 


NC background depends on spectrum 




available; 


and comparable to instrinsic background; 




Signal energy resolution ~ 10%; 


Low Vg signal efficiency (15-20%); 




Underground location 






makes it a multi-purpose detector; 






Cosmic ray rate at 5000ft is ~0. 1 Hz. 






Excellent sensitivity to p ^ 7C^e~^ 


Low efficiency to p ^ K^v 


Liquid 


Technology demonstrated by 


Scale up factor of ~300 is needed; 


Argon 


ICARUS (0.3kT); 




TPC 




Needs considerable R&D for costing; 




Promises high efficiency and 


Not yet demonstrated by 




background rejection; 


simulation of a large detector; 




Has potential to operate 


Needs detailed safety design for 




on (or near) surface; 


deep location in a cavern; 




Could be placed on surface 


Needs detailed demonstration 




either at NuMI Offaxis or DUSEL; 


of cosmic ray rejection; 






Surface cosmic rate ~500kHz; 




Better sensitivity to 


Surface operation limits 






physics program; 



TABLE 11: Comparison of Water Cherenkov to Liquid Argon detector technologies 
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2. INTRODUCTION 



This report details the activities and the results of a several month long study on long baseline 
neutrinos. This workshop (named the US joint study on long baseline neutrinos) was sponsored by 
both Fermi National Accelerator Laboratory and Brookhaven National Laboratory. 

Charge: This study grew out of two parallel efforts. An earlier attempt to create a joint 
FNAL/BNL task force on long baseline neutrinos was initiated by the management of these two 
laboratories. Later the need arose to provide input to the neutrino scientific advisory commit- 
tee (NuSAG) which was asked to address the APS study's recommendation for a next generation 
neutrino beam and detector configuration. The NuSAG charge is in Appendix B. The APS study 
report can be obtained from http://www.aps.org/neutrino/. The study principals created a charge 
with specific scenarios for an accelerator based program. The charge from the chairs of the study 
is in Appendix C. 

Membership: Although the study group was asked to mainly focus on a next generation pro- 
gram within the US, participation from the world wide community of particle physicists was 
sought. In particular, physicists engaged in the European equivalent of this study (the Interna- 
tional Scoping Study: http://www.hep.ph.ic.ac.uk/iss/) were kept abreast of our progress. The Ust 
of physicists who participated in this study by either contributing written material, presentations, 
or discussion is at http://nwg.phy.bnl.gov/fnal-bnl/. 

The membership was divided into several subgroups. The accelerator subgroup studied and 
summarized the proton intensities available mainly from FNAL. The neutrino beam subgroup 
summarized the neutrino beam intensities and event rates for various possibilities. The water 
Cherenkov subgroup summarized the current understanding of the conceptual design of such a 
detector as well as the state of the art in simulating and reconstructing events in such a detector. 
The liquid argon detector subgroup studied the capabilities of such a detector as well as the feasi- 
bility of building a detector large enough to collect sufficient numbers of events. The results from 
each of these groups is either in presentations, technical documents prepared in the near past, or in 
technical documents prepared specifically for this study. 

Scope of the work: As specified in Appendix C, the scope of our work was limited to con- 
ventional horn focused accelerator neutrino beams from US accelerator laboratories. It was asked 
that we study a next generation program by placing massive detectors either off-axis on the sur- 
face for the NuMI beam-line at FNAL, or by building a new intense beam-line aimed towards a 
new deep underground science laboratory (DUSEL) in the western US. The detector technology to 
be considered was either a water Cherenkov detector or a liquid Argon time projection chamber. 
The international scoping study (ISS) on the other hand focused on new technology ideas such as 
beta-beams and muon storage ring based neutrino factories. 
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• In the following we will refer to the NOvA program using the NuMI off axis beam as 
Phase-I. We will not study or comment on this phase extensively since it has been previously 
reviewed extensively, but it will be necessary for us to use the extensive existing material for 
this phase to study the next two items. 

• An upgraded off-axis program with multiple detectors, including a massive liquid argon 
detector, as Phase-II(option A). There could be various versions of Phase-II(option A), with 
or without a liquid argon detector, with a water Cherenkov detector, and/or detectors at 
various locations off axis. We will attempt to elaborate on all of these. 

• A program using a new beam-line towards DUSEL, housing a massive multipurpose detec- 
tor, either a water Cherenkov or a liquid argon detector, will be called Phase-II(option B). 
We will provide information on the DUSEL candidate sites as well as the two options for a 
multipurpose detector. 

Schedule: The study followed the schedule outlined in Appendix |Fj The first meeting of the 
FNAL and BNL management that led to the study was held at BNL on November 14, 2005. The 
charge of the workshop which defined the scope of the work was finalized after the meeting on 
March 6-7, 2006. It was decided at this meeting that since the time for the report was short, it was 
best to create small subgroups to work on individual papers for the study. These papers would be 
distributed to the study group as well as the NuSAG committee as they were prepared. 

A set of presentations were made to the NuSAG committee on May 20, 2006. Results from on- 
going work was reviewed at this meeting. We selected July 15, 2006 as a deadline for preparation 
of the individual papers. Many, but not all, papers were prepared by July 15, and were distributed 
by web-site (http://nwg.phy.bnl.gov/fnal-bnl). 

After discussion within the working group a summary report (this report) was commissioned. 
The contents of this report were reviewed by the study group on September 16-17, 2006. The 
deadline for delivering a preliminary report to NuSAG was October, 2006. 

3. PHYSICS GOALS OF A PHASE-II PROGRAM 

There is now an abundance of evidence that neutrinos oscillate among the three known flavors 
Ve, and v^, thus indicating that they have masses and mix with one another^. Indeed, mod- 
ulo an anomaly in the LSND experiment, all observed neutrino oscillation phenomena are well 
described by the 3 generation mixing 
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with |v, >, i = 1,2,3, the neutrino mass eigenstates. 

Atmospheric neutrino oscillations are governed by a mass squared difference Am^2 = ni^ — ty^ — 
±2.5 X 10^^eV^[5] and mixing angle 023 — 45°; findings that have been confirmed by accelerator 
generated neutrino beam studies at Super- Kamiokande and MINOS [i6ii7i|. 

As yet, the sign of Amgj is undetermined. The so-called normal mass hierarchy, nij, > mj, 
suggests a positive sign which is also preferred by theoretical models. However, a negative value 
(or inverted hierarchy) can certainly be accommodated, and if that is the case, the predicted rates 
for neutrino-less double beta decay will likely be larger and more easily accessible experimentally. 
Resolving the sign of the mass hierarchy is an extremely important issue. In addition, the fact that 
023 is large and near maximal is also significant for model building. Measuring that parameter with 
precision is highly desirable. 



In the case of solar and reactor neutrino oscillations [jSl |9l [TOl, one finds Amjj = nv^ 



8 X lO^^eV^ and 612 — 32°. Again, the mixing angle is relatively large (relative to the analogous 
Cabbibo angle ~ 13° of the quark sector). In addition, Amjj is large enough, compared, to Am32, to 
make long baseline neutrino oscillation searches for CP violation feasible and could yield positive 
results, i.e. the stage is set for a future major discovery (CP violation in the lepton sector). 

Currently, we know nothing about the value of the CP violating phase 5 (0 < 5 < 360°) and 
only have an upper bound ifTTI on the as yet unknown mixing angle 613 (613 < 13°) 



sin^20i3 <0.2 

The value of 613 is likely to be determined by the coming generation of reactor Vg disappearance 
and accelerator based v^u Vg appearance experiments if sin 2 613 > 0.01. Knowledge of 613 
and 8 would complete our determination of the 3 generation lepton mixing matrix and provide a 
measure of leptonic CP violation via the Jarlskog invariant. 
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Jcp = -Sm20i2sin20i3sm2023cos0i3sin5. 
8 

If we use the above limit for 613 then j^p^""^^ < o.05 x sin 5, which could easily turn out to be 
much larger than the analogous quark degree of CP violation J^p"'^'^'' c:=^3 x 10^^. 

Based on our current knowledge and future goals, a phase II neutrino program should include: 

• Completing the measurement of the leptonic mixing matrix, 

• Study of CP violation, 

• Determining the values of all parameters with high precision including Jcp as well as the 
sign of Am32 , 

• Searching for exotic effects perhaps due to sterile neutrino mixing, extra dimensions, dark 
energy etc. 

Of the above future neutrino physics goals, the search for and study of CP violation is of primary 
importance and should be our main objective for several reasons which we briefly outline. 

CP violation has so far only been observed in the quark sector of the Standard Model. Its dis- 
covery in the leptonic sector should shed additional light on the role of CP violation in Nature. Is 
it merely an arbitrary consequence of inevitable phases in mixing matrices or something deeper? 
Perhaps, most important, unveiling leptonic CP violation is particularly compelling because of its 
potential connection with the observed matter-antimatter asymmetry of our Universe, a funda- 
mental problem at the heart of our existence. The leading explanation is currently a leptogenesis 
scenario in which decays of very heavy right-hand neutrinos created in the early universe give rise 
to a lepton number asymmetry which later becomes a baryon-antibaryon asymmetry via the B-L 
conserving 't Hooft mechanism of the Standard Model at weak scale temperatures. 

Leptogenesis offers an elegant, natural explanation for the matter-antimatter asymmetry; but 
it requires some experimental confirmation of its various components before it can be accepted. 
Those include the existence of very heavy right-handed neutrinos as well as lepton number and 
CP violation in their decays. 

Direct detection of those phenomena is highly unlikely; however, indirect connections may be 
established by studying lepton number violation in neutrinoless double beta decay and CP violation 
in ordinary neutrino oscillations. Indeed, such discoveries will go far in establishing leptogenesis 
as a credible, even likely scenario. For that reason, neutrinoless double beta decay and leptonic CP 
violation in neutrino oscillations are given very high priorities by the particle and nuclear physics 
communities. 
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Designing for CP violation studies in next generation neutrino programs has other important 
benefits. First , the degree of difficulty to establish CP violation and determine f^P*°"^^^ is demand- 
ing but doable. It requires an intense proton beam of about 1-2 MW and a very large detector 
(250 ~ 500 kton Water Cherenkov or a Uquid argon detector of size ~ 100 kTon which could be 
equivalent in sensitivity due to its better performance). Such an ambitious infrastructure will allow 
very precise measurements of all neutrino oscillation parameters as well as the sign of Am32 via 
V;u disappearance and v^i Vg appearance studies. It will also provide a sensitive probe of 
"New Physics" deviations from 3 generation oscillations, perhaps due to sterile neutrinos, extra 
dimensions, dark energy or other exotic effects. 

A well instrumented very large detector, in addition to its accelerator based neutrino program, 
could be sensitive to proton decay which is one of the top priorities in fundamental science. As- 
suming that it is located underground and shielded from cosmic rays, it can push the limits on 
proton decay into modes such as p — > e^n^ to lO^^yr sensitivity or beyond, a level suggested by 
gauge boson mediated proton decay in super-symmetric GUTs. Indeed, there is such a natural 
marriage between the requirements to discover leptonic CP violation and see proton decay (i.e. an 
approximately 500 kTon water Cherenkov detector) that it could be hard to imagine undertaking 
either effort without being able to do the other. 

Such a large detector would also have additional physics capabilities. It could study atmospheric 
neutrino oscillations with very high statistics and look for the predicted relic supernova neutrinos 
left over from earlier epochs in the history of the Universe, a potential source of cosmological 
information. Also, if a supernova should occur in our galaxy (expected about every 30 years), such 
a detector would see about 100,000 neutrino events. In addition, it could be used to look for signals 
oin — h oscillations in nuclei and highly penetrating GUT magnetic monopoles which would leave 
behind a trail of monopole catalyzed proton decays. 

The physics potential of a very large underground detector is extremely rich. The fact that it 
can also be used to determine (or bound) leptonic CP violation and measure all facets of neutrino 
oscillations gives such a facility outstanding discovery potential. It would be an exciting, central 
component of the world's particle physics program for many decades. On the other hand, a staged 
approach using existing beam facilities should also be explored to determine an optimum strategy. 

4. STRATEGIES FOR THE PHASE-II PROGRAM USING A CONVENTIONAL BEAM 

In this section we will describe the essential features of an off-axis narrow band beam versus an 
on-axis broad band beam. We will then briefly summarize how these features can be used to extract 
the CP violation effect as well as all the other parameters of importance in neutrino oscillations. 

Throughout this report we are concerned with conventional horn focused beams in the US: 
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the existing NuMI beam at FNAL or a new super neutrino beam that could be optimized for a 
detector at a new deep underground national laboratory (DUSEL) with a possible large detector 
(either underground or on the surface). The measurement of most interest is always the appearance 
measurement, Vg, for which the horn focused beam has a limitation from the irreducible 

background of Vg contamination in the beam. The level of contamination depends on neutrino 
energy and also the beam design and the off-axis angle, but it is in the range of ~ 0.5 — 1% for 
most practical beams. This contamination comes from decays of muons and kaons in the beam. 
These cannot be completely eliminated. The second source of background is neutral current events 
that mimic electron showers. This background is considered reducible by detector design. In 
particular, a fine grained detector such as a liquid argon TPC detector will be capable of reducing 
such background to very small levels. Most of the remaining report will be concerned with the best 
strategy for obtaining sufficient signal events while reducing these backgrounds. In this section we 
will not discuss the issues of backgrounds in detail, but give a guide to the signal spectra, event 
rates and comment on the implications. 

Figures [T] [2} [3| and |4] show the spectra of concern. Care is required in comparing these plots 
because they are plotted on a logarithmic energy scale. The normalization is per GeV of neutrino 
energy per kTon of detector mass per MW x 10^ sec protons of the appropriate energy on target. 

These spectra were obtained by detailed simulations using the GNuMI computer program[7J. 
For these figures a simple recipe was used to obtain charged current event rate [[Till : a cross section 
of 0.8 X 10^^^ ctrp-/ GeV (0.35 x lO^^^cm^/GeV for anti-neutrinos) was used above 0.5 GeV and 
the quasi-elastic cross section was used below 0.5 GeV. There could be small differences due to 
the detector target type (water, argon, etc.), but this is a good approximation [fT3l . For figures [T] 
and |2] we have used the low energy (LE) setting of the NuMI beam configuration which gives a 
better flux at the 40 km site. Reference [[T2| contains spectra for other choices. For all the off-axis 
spectra 120 GeV protons were used and the normalization is for MW x 10^ sec protons; for 120 
GeV protons this corresponds to 5.2 x 10^*^ protons. 

For Figures |3] and |4[ the GNuMI program was modified for a wide band low energy (WBLE) 
design for the horns as well as a new decay tunnel with 4 m diameter and 400 m length; these 
are described in detail in [fT4l . For the WBLE beam, there is a choice of running with protons 
from 40 GeV to 120 GeV. For these plots we have chosen 60 GeV protons. The normalization is 
for MW X lO^^ec protons of 60 GeV. The spectra shown here should not be considered optimum. 
After thorough design and optimization there could be modest improvements, but at this point we 
are confident that these numbers are sufficiently good for this review. 

For Figures[T|to|4]we have superimposed the expected probability of — > Vg conversion for the 
appropriate distance and for the following oscillation parameters: Am^2 = 0.0025eV^, Am^^ = 8 x 
10-^ eV^, sin^ 2012 = 0.86, sin2 2023 = 1-0, and sin2 20i3 = 0.04; the curves are for several choices 
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of the CP phase and the left and right hand side plots are for the two different mass orderings. 



In Table III we have calculated the rate of electron appearance events for various scenarios by 
integrating the spectrum together with the appearance probability. This event rate is for all charged 
current events; no detector efficiency factors are applied. A detector with efficient fiducial mass of 
100 kTon is assumed with ~ 10^ sec of running time with 1 MW of proton beam. No consideration 
for backgrounds, energy thresholds, or resolution effects are in this table. 

Also note that Figures [T] to |4] do not show the event rates from anti-neutrinos. These can be 
obtained from the study web-site [ITSl . We include anti-neutrino rates and spectra in later sections 
with more detail. We have included anti-neutrino event rate in Table Hill 

After considering the figures and the table we make the following observations: 

• For simplicity we look at the electron neutrino event rate at dcp — and compare it to dcp = 
—90". In the limit that one has resolved the mass hierarchy using the anti-neutrino data, the 
modulation of the neutrino rate with 5cp will give us the CP parameter measurement that we 
seek. One can immediately see that the size of the CP effect for the maximum CP (— 90*^) 
is approximately 3c7. To achieve this within a year of running (with no consideration for 
efficiencies, backgrounds, etc.) the efficient fiducial mass of the detector must be lOOkTon 
range if the accelerator power is limited to be ~1 MW. This conclusion is regardless of the 
eventual choice for the beam-line. 

• The size of the CP effect (for the maximum 90") increases modestly from ~ 3a for the off 
axis (810 km) 12 km option to about 4.8a for 2500 km. Much of this increase can be traced 
to the large CP effect at the higher oscillation nodes that become available for the larger 
distances. The loss of statistics due to distance (as l/L^) is largely compensated by the 
increase in the strength of the CP related signal [fT6l [TTl [IHl. By combining the 40 km off- 
axis rates with the 12 km there is also a modest improvement in the overall CP measurement. 
Nevertheless, for the choice of spectra in this report, the baseline length related effects for a 
CP measurement are not dramatic for the range of choices in this study. 

• Remarkably, it should also be noticed that the size of the CP effect in the number of sigma is 
approximately the same for the different values of sin 20i3. It has been pointed out, there- 
fore, that for sin^20i3 > 0.003, which is the range accessible for conventional accelerator 
beams, the size of the exposure (efficient fiducial mass multiplied by the total incident beam 
power) needed to obtain a good measurement of the CP parameter is independent of 

\. This is explained by the following argument. The asymmetry defined by 
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is proportional to J^p and therefore grows linearly with sin 613, but /'(v^ Vg) is to 
leading order proportional to sin 20i3 and therefore the statistical figure of merit, the error 
on the asymmetry A should have little dependence on 613. 

The size of the matter effect (the difference between the event rate for the two choices of 
mass ordering) is approximately 3(7 for the 12 km off axis location for sin^20i3 = 0.02 for 
neutrino rates alone. It is a much larger effect for longer baselines. The probability curves 
show that the effect is large for the first oscillation node in all cases. This effect will clearly 
compete with the CP effect and must be determined along with the CP effect for clarity. The 
matter effect clearly is much stronger for larger value of 613, and therefore for a larger value 
of 013, it will be easier to determine the mass hierarchy. 

Examination of the probability curves in Figure [T] shows that the 12 km off axis spectrum is 
sensitive mainly to the first oscillation node. The probability is affected not only by the CP 
phase, but also by the value of 813, the mass ordering, the uncertain values of other param- 
eters such as Am^2 ^23- Also note that the probability curves at any particular energy 
have degeneracies in the CP phase. These degeneracies have been discussed in the literature 
[|2Tll22ll23l . Therefore, to make a clean determination of CP violation, one either needs very 
good energy resolution (to exploit the small energy dependence within the first node) with 
good statistics, or one needs to perform another measurement at the high oscillation node by 
placing another detector further off- axis. This is one of the options to be examined in this 
report. 

Examination of the probability curves in Figures |3] and |4] shows that the energy dependence 
of the probability can be measured in a single detector by creating a beam spectrum that 
matches the first few nodes over the > 1000 km long baseline. Obviously, in such a sce- 
nario the neutrino energy must be measured in the detector with sufficient resolution while 
suppressing backgrounds [flTll . This is also an option to be considered in this report. An 
illustration of how the various degeneracies affect the measurement is shown in Figure |5] 
The figure illustrates the energy dependence for neutrino running only. It is clear that nar- 
row band running will have additional ambiguities. How these can be broken with additional 
anti-neutrino running or with high statistics and resolution will be discussed later. 

The neutrino event rate is roughly proportional to the total proton beam power; the exact 
numbers and deviations from this rule will be discussed below. The total power that can be 
obtained from FNAL Main Injector after upgrades increases with the output proton energy, 
and therefore it is important to maintain the highest possible proton energy for either the 
off-axis or on-axis scenarios. For the off-axis experiment the preferred running is at the 
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FIG. 1: (in color) Spectrum of charged current events at a 12 km off-axis location at 810 km on the 
NuMI beam-line. The spectrum is normalized per GeV per MW x 10^ sec protons of 120 GeV. The low 
energy (LE) setting of the NuMI beam-line is used for this plot. Overlayed is the probability of Ve 
conversion for sin^ 20i3 = 0.04 with rest of the oscillation parameters as described in the text. The left plot 
is for regular mass ordering and right hand side is for reversed mass ordering. Figure includes no detector 
effects such as efficiencies, resolution, or backgrounds. 

highest, 120 GeV, proton energy. For the FNAL-to-DUSEL option, there could be significant 
advantage at running with lower proton energy. This will reduce the long high energy tail 
> 5GeV of the neutrino spectrum. This tail is outside the interesting oscillation region and 
may contribute increased background in the form of neutral current events that reconstruct 



to have lower neutrino energy. The event rates given in table III for WBLE assume running 
with 1 MW of power at 60 GeV. In the following we will comment on how 1 MW power can 
be obtain while maintaining the a flux with low high energy neutrino tail. The easiest way, 
of course, is by having a small off-axis angle. The flux that could be obtained with a 0.5" 
off-axis angle to DUSEL at 1300 km is shown in Figure [6) 

We will now explore the above observations in further detail including the feasibility of beams 
and detectors, current best knowledge on the performance of detectors, and requirements for other 
physics related applications of these very large detector facilities. 
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FIG. 2: (color) Spectrum of charged current events at a 40 km off-axis location at 810 km on the NuMI 
beam-line. The spectrum is normalized per GeV per MW x 10^ sec protons of 120 GeV. The low energy 
(LE) setting of the NuMI beam-line is used for this plot. Overlayed is the probability of — > Vg conversion 
for sin^ 2013 = 04 with rest of the oscillation parameters as described in the text. The left plot is for regular 
mass ordering and right hand side is for reversed mass ordering. Figure includes no detector effects such as 
efficiencies, resolution, or backgrounds. 



5. ACCELERATOR REQUIREMENTS 



All phases of the envisioned US neutrino accelerator program, Phase-I(NOvA), Phase-II(option 
A), or Phase-II(option B), require upgrades to the existing proton accelerator infrastructure in the 
US. Phase-I upgrades, already planned at FNAL, will increase the Main Injector extracted beam 
power to 0.7 MW at 120 GeV (this is called "proton plan-2" and has been incorporated in the 
NOvA project). The plan to further upgrade the Main Injector to 1.2 MW is called "the SNuMI 
Project" [[24ll . Phase-II will benefit from these upgrades. 

We have used beam power in the range of ~0.5 to 2 MW for high energy protons (>30 GeV) 
in our calculations because this level of beam power is now considered the next frontier for current 
accelerator technology [l25l |26l |271 |28l |29l |30l and also necessary to obtain sufficient event rate 
to perform the next stage of neutrino oscillation physics. The technical limitations arise from the 
need to control radiation losses, limit the radiation exposure of ground water and other materials, 
and the feasibility of constructing a target and horn system that can survive the mechanical and 
radiation damage due to high intensity proton pulses [29]. 

We quote event rates either in units of MW x IQp sec or number of protons on target (POT). A 
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FIG. 3: (color) Spectrum of charged current events using a new wide band beam from FNAL to a location 
at 1300 km. The spectrum is normalized per GeV per MW x 10^ sec protons of 60 GeV. Overlayed is the 
probability of — > Vg conversion for sin^20i3 = 0.04 with rest of the oscillation parameters as described in 
the text. The left plot is for regular mass ordering and right hand side is for reversed mass ordering. Figure 
includes no detector effects such as efficiencies, resolution, or backgrounds. 

convenient formula for conversion is below. 

2Q. 1000 X BeamPower{MW) x r(10^^) 



POT(\0^ 



1.602 xEp{GeV) 

where T is the amount of exposure time in units of 10^ s and Ep is the proton energy. We now 
briefly summarize the understanding of high energy proton beam power at the two US accelerator 
laboratories where high intensity proton synchrotrons are operational, Fermilab and Brookhaven. 

FNAL Main injector (MI): Discussion is currently underway to increase the total power from 
the 120 GeV Main Injector (MI) complex after the Tevatron program ends [|24ll25l . In this scheme 
protons from the 8 GeV booster, operating at 15 Hz, will be stored in the recycler (which becomes 
available after the shutdown of the Tevatron program) while the MI completes its acceleration 
cycle, which is shortened from the current 2.2 sec to 1.33 sec. In a further upgrade the techniques 
of momentum stacking using the antiproton accumulator, and slip-stacking using the recycler will 
raise the total intensity in the MI to ~ 1 .2 MW at 120 GeV [|26l|. In the rest of this report this will be 
called the SNuMI plan. In the ideal case, the length of the acceleration cycle is proportional to the 
proton energy, making the average beam power proportional to the final proton energy. However, 
fixed time intervals in the beginning and the end of the acceleration cycle are required for stable 
operation. These become important at low energies and reduce the performance below the ideal. 
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FIG. 4: (color) Spectrum of charged current events using a new wide band beam from FNAL to location 
at 2500 km. The spectrum is normalized per GeV per MW x 10^ sec protons of 60 GeV. Overlayed is the 
probability of — > Vg conversion for sin^20i3 = 0.04 with rest of the oscillation parameters as described in 
the text. The left plot is for regular mass ordering and right hand side is for reversed mass ordering. Figure 
includes no detector effects such as efficiencies, resolution, or backgrounds. 



Current projections suggest that ~ 0.5 MW operation between 40 — 60 GeV and > 1 MW operation 
at 120 GeV is possible. 

More ambitious plans at FNAL call for replacing the 8 GeV booster with a new super- 
conducting LINAC that can provide 1.5 x 10^^ H ions at 10 Hz corresponding to 2 MW of 
total beam power [|27]| . Some of the 8 GeV ions could be injected into the MI to provide high 
proton beam power at any energy between 30 and 120 GeV; e.g., 40 GeV at ~ 2 Hz or 120 GeV at 
~ 0.67 Hz. Such a plan allows for flexibility in the choice of proton energy for neutrino production. 
This plan will be called the high intensity neutrino source upgrade (HINS). 

The projected proton intensity from the main injector for the successive upgrades at FNAL is 
shown in Figure [7 J|28]| . A reviewed cost estimate that has been included in the NOvA project for 
the 700 MW (proton plan-2) upgrade is $60M. The cost of the complete SNuMI plan (to 1.2 MW) 
is at the moment very preliminary at ~ $54M (without overhead or contingency factors). The 
HINS upgrade is estimated to be approximately >$300M. 

BNL AGS: The BNL Alternating Gradient Synchrotron (AGS) operating at 28 GeV currently 
can provide about 1/6 MW of beam power. This corresponds to an intensity of about 7 x 10^-^ 
protons in a 2.5 microsecond pulse every 2 seconds. The AGS complex can be upgraded to provide 
a total proton beam power of 1 MW ll30l . The main components of the accelerator upgrade at BNL 
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FIG. 5: (in color) Spectrum of charged current Ve events using the LE beam tune at 12 km off-axis 810 
km location (left) and with a new wide band beam from FNAL (using 60 GeV protons) to a location at 1300 
km. The spectra are normalized for 600MW x 10^ ^ec and the width of the band indicates the statistical error. 
The parameters used for oscillations are shown in the figure, the remaining parameters are as described in 
the text. Figure includes no detector effects such as efficiencies, resolution, or backgrounds. 




FIG. 6: (in color) Spectrum of charged current events using a new wide band beam from FNAL to 
location at 1300 km with slightly off axis location (12km) to reduce the high energy tail. The spectrum 
is normalized per GeV per MW x lOp sec protons of 120 GeV. Overlay ed is the probability of Vg 
conversion for sin^20i3 = 0.04 with rest of the oscillation parameters as described in the text. The left plot 
is for regular mass ordering and right hand side is for reversed mass ordering. Figure includes no detector 
effects such as efficiencies, resolution, or backgrounds. 
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TABLE ni: This table contains signal event rates after Ve (also for anti-neutrinos) conversion for the 
various scenarios described. The event rates here have no detector model or backgrounds. The units are 
charged current events per 100 kTon of detector mass for 1 MW of beam for 10^ sec of operation. For NuMI 
running we assume 120 GeV protons in the LE tune and for WBLE we have assumed 60 GeV protons. The 
charged current cross sections applied as well as the oscillation parameters used are described in the text. 



are a new 1.2 GeV super-conducting LINAC to provide protons to the existing AGS, and new 
magnet power supplies to increase the ramp rate of the AGS magnetic field from about 0.5 Hz to 
2.5 Hz. For 1 MW operation 28 GeV protons from the accelerator will be delivered in pulses of 
9 X 10^^ protons at 2.5 Hz. It has been determined that 2 MW operation of the AGS is also possible 
by further upgrading the synchrotron to 5 Hz repetition rate and with further modifications to the 
LINAC and the RF systems. The AGS 1 MW upgrade is estimated to cost $343M (TEC) including 
contingency and overhead costs. This cost has been reviewed internally at BNL. 
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FIG. 7: Proton beam power from the Fermilab main injector as a function proton energy for various scenar- 
ios. Lowest (blue) curve is for the current complex running concurrently with the Tevatron. Second (green) 
curve is for the proton plan upgrades, third (light blue) curve is for SNuMI recycler stage which will take 
place after the termination of the Tevatron program, fourth (red) is for the accumulator stage upgrades, the 
uppermost (brown) is for the HINS upgrade which calls for a new 8 GeV LINAC injector. 

6. TARGET AND HORN DEVELOPMENT 

All phases of the envisioned US neutrino accelerator program, Phase-I (NOvA), Phase- 
II(option A), or Phase-II(option B), require substantial development for a new target capable of 
operating at high proton intensities and perhaps new focusing horn optics. 

Current understanding of targets, and R&D in progress is summarized in [|29l[3T1l . The neutrino 
event rate is approximately proportional to the total proton beam power (energy times current) 
incident on the target. The parameters for target design to be considered for a given power level 
are proton energy, pulse duration, and repetition rate. In addition to these the shape and size of the 
beam spot on the target, and the angle of incidence could also be varied. Studies over the last few 
years have come to the acceptance that with optimal choice of the above variables the upper limit 
for a solid target operation is ~ 2 MW. For a given accelerator facility these parameters tend to be 
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correlated and constrained, and therefore a practical limit for a solid target with current technology 
is probably between 1 and 2 MW. Nevertheless, considerable work is needed to achieve a practical 
design for such a high power solid target and its integration into a focusing horn system. Above 
2 MW, liquid targets are likely the better choice, but these devices will require considerable R&D 
and testing before they can be considered practical. 

Target R&D which includes understanding of materials as well as engineering issues of inte- 
gration is a critical item for the physics program considered in this report. 

7. NEUTRINO BEAM-LINES 

There is currently good experience in building and operating high intensity neutrino beam-lines 
in the US. The study group has concluded that it is possible to use an existing or build a new super 
neutrino beam-line based on current technology or extensions of current technology and operate it 
for the physics program described in this report. 

In the following we summarize the status of US high energy accelerator neutrino beam-lines. 
There are two additional accelerator neutrino beam-lines in the world with comparable technical 
requirements: the CERN to Gran Sasso neutrino beam which is now operating, and the JR^RC to 
Super- Kamiokande neutrino beam which will start operation in a few years. We will not report on 
these in this report, but this study included presentations from these facilities. It is clear that there 
is plenty of communication and shared technical information between these centers and the US. 

7.1. NuMI 

The design and operation of the NuMI beam-line was reported in [31]. In the NuMI beam 120 
GeV protons from the Main Injector, in a single turn extraction of ~ lOjis duration every ~ 2sec, 
are targeted onto a 94cm long graphite target. A conventional 2 horn system is used to charge 
select and focus the meson beam into a 675 m long, 2 meter diameter evacuated decay tunnel. The 
NuMI beam-line is built starting at a depth of ~ 50m and is aimed at a downwards angle of 3.3 
deg towards the MINOS detector in Minnesota at a distance of 735 km from the production target. 
The flux of the resulting neutrino beam is well known and will be described in a separate section 
below. 

NuMI beam transport, target, horns, and shielding were designed for operation with 4 x 
10^^ protons / pulse with a beam power of 400 kW. The goal is to average 3.7 x \0^^ protons / year . 
The first year run of NuMI achieved typical beam intensities of 2.5 x \0^^ protons/ pulse or 200 
kW. The total integrated exposure was 1 .4 x 10^*^ protons on target for the period from March 2005 
to March 2006. A number of technical problems were encountered and solved during this time: at 
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the start of the run the cooling water hne to the target failed, one of the horns had a ground fault, 
and most notably a detailed study of the tritium production from the beam-hne had to be carried 
out. Various monitoring systems as well systems to collect tritiated water were installed to elimi- 
nate the amount of tritium going into cooling water and the environment. The experience gained 
from NuMI operations is indeed invaluable for future operation of neutrino beams. 

NuMI beam- line was built at a total cost of $109M (TEC). The construction time was approx- 
imately ~5 years. The beam- line became operational in March of 2005. Upgrade and operation 
of the NuMI beam-line for higher intensity for Phase-I is included in the new SNuMI conceptual 
design report at FNAL [|24l |28l . It is anticipated that for operation at 1 MW, the primary proton 
beam- line, the target and horns, and cooling systems in the target hall will require upgrades. New 
He bags and upgrades to the high radiation work areas will also be installed. The total preliminary 
cost of this upgrade (~ lOM) is included in the cost of the proton plan-2 upgrade described in 
Section m 

7.2. Beam towards DUSEL 

Members of this study group [l28l have examined the possible siting and construction of a new 
beam-line towards one of the site candidates for DUSEL, either Henderson mine in Colorado or 
Homestake mine in South Dakota. The study group has concluded that there are no technical 
limitations to building such a beam-hne on the Fermilab site using the same extraction line from 
the main injector as the NuMI beam-line. The study group has found significant advantage for 
lower energy neutrino flux in making the diameter of the decay tunnel for the new beam-line up to 
4 meters. 

The new beam-line at FNAL would use the same extraction from the Main Injector into the 
NuMI line; a new tunnel would pick up the proton beam from the present tunnel and transport it in 
the western direction with the same radius of curvature as the Main Injector so that up to 120 GeV 
protons can be used with conventional magnets. There is adequate space on the Fermilab site to 
allow a new target hall with 45 m length and a decay tunnel of length 400 meter and diameter of 
4 meters. This will allow the location of a near detector with ~ 300 meters of length from the end 
of the decay pipe. The new decay pipe would point downwards at an angle of 5.84" to Homestake 
(1289km from FNAL) or 6.66'' to Henderson (1495 km from FNAL). The diameter of the decay 
tunnel is a crucial parameter for both the neutrino beam intensity and the cost and feasibihty of 
the beam-line; it will require detailed optimization. With our present understanding, construction 
of a 4 meter diameter decay tunnel with adequate shielding for eventually 2 MW of operation is 
possible. If the additional concrete shielding is found to be inadequate then the decay pipe would 
have to be reduced to 3 meter diameter because of the maximum possible span of excavation in the 
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rock under FNAL ll32l . The thickness of the shielding has been scaled from the NuMI experience, 
but the implications of the wider diameter for radiation issues (in particular, tritium production) 
will need careful study. After optimization, the cost of such a project can be reliably estimated 
from the known cost of the NuMI project. 

A new beam-line from BNL-AGS to either Homestake (2540 km) or Henderson (2770 km) has 
also been examined in a BNL report[30J. They have made the choice of building the beam- line 
on a specially constructed hill where the shielded target station is located on top of the hill and 
the meson decay tunnel is on the downward slope of the hill pointing towards DUSEL at an angle 
of 11.7" (Homestake) or 13.0" (Henderson). Due to the limitations on the height of the hill, the 
decay tunnel length is restricted to be ~200 meters with a diameter of 4 meters. The cost of such a 
beam- line including construction of the hill and proton transport to the top of the hill was estimated 
to be $64M (TEC) including contingency and overhead; this cost has been reviewed internally at 
BNL. Further work on this option has not been part of this study. 

8. EVENT RATE CALCULATIONS 

The neutrino flux and the numbers of expected events with and without oscillations were cal- 
culated for both the NuMI off-axis beam and a new broadband beam towards DUSEL. This calcu- 
lation assumes no detector resolution model or background rejection capabihty. Both calculations 
were performed using the same GEANT based GNuMI code. This code has been extensively 
tested as part of the MINOS collaboration. It has been verified against recent data in the MINOS 
near detector. The code and associated cross section model is known to produce agreement with 
the MINOS near detector event rate per proton to about 10% at the peak of the spectrum and of 
the order of 20-30% in the tails of the spectrum with no adjustments. We have also calculated 
anti-neutrino event rates. The accuracy here is worse simply because of the lack of data from the 
NuMI beam-line. The anti-neutrino spectra have disagreements between various production codes 
of ~ 30%. We believe this is sufficient accuracy for the purposes of this study. 

It is very likely that neither the specific off-axis configuration nor the broad-band configura- 
tion is highly optimized for the physics under consideration. Such optimization could result in 
modest gains, especially at low energies. At this stage there is good confidence that the possible 
improvements will not change the overall picture and sensitivity outlined in this report. 

8.1. NuMI off-axis locations 

We have calculated the neutrino flux and event rates at various off-axis distances from the NuMI 
beam-line. NuMI was assumed to be configured in the medium energy (ME) or low energy (LE) 
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beam configuration for the results quoted here. The low energy configuration provides better event 
rate at the 40 km off-axis location in the low energy peak. There is, however, event rate loss at the 
12 km location. 



The details of the calculation, as well as the spectra are in [[T2J. Tables IV for neutrino running 
and |V] for antineutrino running summarizes these event rates. The normalization is per MW x 
10^ sec protons of 120 GeV and for 1 kTon of efficient detector mass. There are no corrections for 
the type of target nucleus in the detector. There are no efficiencies for reconstruction or fiducial 
cuts in this calculation. 

We have used tabulated cross sections to calculate the event rates in the various columns. The 
column labeled "v^ CC" is the total charged current muon neutrino event rate, "v^ CC osc" is the 
charged current muon neutrino event rate after oscillations. "Vg CC beam" is the charged current 
rate of electron neutrino contamination in the beam. "Vg QE beam" is the charged current quasi- 
elastic event rate of electron neutrino contamination in the beam. "NC-l;r^" is the rate of neutral 
current single pion production integrated over the noted energy range; no detector related rejection 
is assumed in this table, "v^x — > Vg CC" is the charged current event rate of electron neutrinos after 
oscillations using the oscillations parameters described in Section |4| "v^ Vg QE" is the quasi- 
elastic rate of electron neutrinos after oscillations using the oscillations parameters described in 
Section |4j For example, the total CC event rate in 5 years with 1.7 x 10^ sec/yr in a 100 kton 
detector without oscillations at 40 km (LE) off axis can be calculated to be 5.38 x 100 x 5 x 1.7 = 
4573. This event count includes events from both the pion and the kaon peaks at about 0.5 and 4 
GeV, respectively. 



8.2. Wide band beam towards DUSEL 

The spectra and the event rate for a beam towards DUSEL were calculated by using the same 
GNuMI framework but the geometry of the target, horns, and the decay tunnel was changed. The 
full calculation and the resulting spectra are described in [fT4l . The integrated event rates are shown 



in Table VI and Table VII There are a few of comments of importance: 



The calculation in the table is for 1300 km (the FNAL to Homestake distance), but it could 
be easily converted to 1500 km (the distance to Henderson). The unoscillated rate scales 
as 1/r^, but the oscillated event rate scales according to the oscillation function. When we 
demonstrate the full sensitivity calculation later in this report we include the variation with 
distance. For 1300 versus 1500 km this variation is small. 



Ealier work on sensitivity used a 28 GeV proton beam [|20ll . The total CC event rate in 
100 kTon efficient fiducial mass after 5 years at 1.7 x 10^ sec/yr without oscillations using 



29 



TABLE IV: Signal and background interaction rates for various NuMI beam configurations, baselines and 
off-axis distances. Rates are given per MW.lO^s.kT. The rates are integrated over the range 0-20 GeV. For 
— > Ve oscillations a value of sin^20i3 = 0.04 and Amgj = 2.5 x 10^^ eV^ is used. No detector model is 
used. 



Distance off-axis 


V;:CC 


CC osc 


CC beam 


Ve QE beam 


NC-l7r° 


V;, ^ Ve CC 


Ve QE 


NuMI LE tune at 700 km 


0km 


400.2 


267.6 


4.55 


0.444 


21.2 


3.66 


0.676 


40 km 


4.81 


2.66 


0.190 


0.047 


0.525 


0.071 


0.038 


NuMI LE tune at 810 km 


0km 


299.0 


187.4 


3.40 


0.332 


15.8 


3.10 


0.551 


6 km 


198.6 


107.0 


2.59 


0.275 


11.9 


2.53 


0.506 


12 km 


84.4 


31.9 


1.57 


0.193 


6.79 


1.41 


0.367 


30 km 


11.6 


8.38 


0.353 


0.070 


1.32 


0.107 


0.046 


40 km 


5.38 


2.91 


0.195 


0.045 


0.596 


0.084 


0.045 


NuMI ME tune at 810 km 


0km 


949.1 


781.1 


7.14 


0.485 


30.6 


4.71 


0.527 


6 km 


304.9 


191.4 


3.83 


0.313 


14.9 


3.19 


0.491 


12 km 


80.5 


32.0 


1.81 


0.174 


5.74 


1.33 


0.330 


30 km 


8.59 


5.52 


0.321 


0.051 


0.81 


0.094 


0.038 


40 km 


4.14 


2.40 


0.168 


0.032 


0.427 


0.054 


0.022 



Ep = 28 GeV protons with 1 MW running is 44625 events integrated over 1-20 GeV. It 
should be kept in mind, however, that according to [|25l . the available beam power is less 
for lower energies (see Fig. [7]). In the technical note [fT4ll it has been shown that the 40-60 
GeV spectrum could be very similar to the 28 GeV with considerable increase in event rate 
per unit beam power. It has also been shown that it is possible to run at the full energy of 
120 GeV and still obtain essentially the same spectrum as the 28 GeV one with a small 0.5^^ 
off-axis angle. With such a choice the neutrino (antineutrino) event rate is 76415 (28475) 
for 100 kTon and 5 yrs for 1 MW and 1.7 x 10^ sec/yr. 



Tables VI and VII represent our present understanding of creating such a beam. When 
optimization is performed coupled to the complete understanding detector performance ver- 
sus energy, the spectrum could be adjusted to give the best signal/background performance. 
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TABLE V: Signal and background interaction rates for various NuMI anti-neutrino beam configurations, 
baselines and off-axis distances. Rates are given per MW.lO^s.kT. The rates are integrated over the range 
0-20 GeV. For v^^ Vg oscillations a value of sin^20i3 = 0.04 and Amj^^ = 2.5 x 10"^ eV^ is used. No 
detector model is used. 



Distance off-axis 


v^CC 


CC osc 


Ve CC beam 


Ve QE beam 


NC-l7r° 


^ Ve CC 


Ve QE 


NuMI LE tune at 700 km 


0km 


157.6 


102.3 


1.69 


0.306 


19.3 


1.25 


0.306 


40 km 


1.64 


0.905 


0.063 


0.021 


0.544 


0.024 


0.016 


NuMI LE tune at 810 km 


0km 


117.7 


71.0 


1.26 


0.229 


14.4 


1.026 


0.285 


6 km 


77.6 


39.8 


0.925 


0.179 


10.8 


0.800 


0.241 


12 km 


31.7 


10.9 


0.545 


0.116 


6.29 


0.388 


0.145 


30 km 


3.87 


2.69 


0.122 


0.035 


1.31 


0.043 


0.025 


40 km 


I.8I 


0.97 


0.066 


0.021 


0.609 


0.029 


0.018 


NuMI ME tune at 810 km 


0km 


350.6 


285.1 


2.53 


0.349 


23.6 


1.59 


0.316 


6 km 


112.8 


69.0 


1.28 


0.208 


11.9 


1.011 


0.259 


12 km 


27.7 


9.83 


0.601 


0.105 


4.76 


0.348 


0.125 


30 km 


2.66 


1.67 


0.109 


0.027 


0.70 


0.027 


0.014 


40 km 


1.27 


0.73 


0.057 


0.016 


0.376 


0.015 


0.008 



This could be accomplished by optimizing the horn optics and/or inserting secondary targets 
(plugs) that remove high energy pions from the beams (see ||33l ). 

• We have integrated the rates of various types of events over the same energy interval 0-20 
GeV for Tables |IV] to IVIII It should be understood that there is considerable variation in 
the signal to background ratio as a function of energy. To get a full appreciation of this we 
recommend the reader to explore the spectra at the study web-site [fTSl . The variation also 
depends on oscillation parameters. In particular, it should be noted that the CP violating 
phase as well as the mass hierarchy is responsible for moving the peak of the oscillation 
probability by as much as ~0.5 (0.7) GeV for the 810 (1300) km baseline. This variation 
coupled to the width of the useful spectrum and the detector energy resolution has an impact 
on the parameter sensitivity of the program. 
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TABLE VI: Signal and background interaction rates at 1300 Km (Fermilab-HOMESTAKE) using different 
WBLE beam energies and off-axis angles. The rates integrated over the neutrino energy range of - 20 
GeV. Rates are given per MW.lO^s.kT. For Ve oscillations a value of sin^20i3 = 0.04 and Arn^i = 

2.5 X 10^^ eV^ is used. No detector model is used. 



Degrees off-axis 


V, CC 


CC osc 


v^, CC beam 


v^, QE beam 


NC-Itt" 


V, - V, CC 


— V, QE 


WBLE 120 GeV at 1300 km with decay pipe 2m radius 380 m length 


0° 


198.2 


104.9 


1.89 


0.179 


9.11 


2.85 


0.408 


0.5° 


89.9 


37.9 


1.22 


0.140 


5.62 


1.62 


0.300 


1.0° 


34.2 


19.5 


0.621 


0.095 


2.95 


0.470 


0.129 


2.5° 


4.66 


2.36 


0.116 


0.032 


0.550 


0.094 


0.049 


WBLE 60 GeV at 1300 km with decay pipe 2m radius 380 m length 


0° 


151.0 


69.2 


1.34 


0.169 


7.83 


2.53 


0.403 


0.5° 


77.2 


28.7 


0.906 


0.134 


5.33 


1.52 


0.305 


1.0° 


33.3 


18.4 


0.520 


0.098 


3.08 


0.480 


0.141 


2.5° 


5.05 


2.56 


0.120 


0.035 


0.611 


0.105 


0.058 


WBLE 40 GeV at 1300 km with decay pipe 2m radius 380 m length 


0° 


110.4 


44.4 


1.02 


0.159 


6.50 


2.05 


0.357 


WBLE 28 GeV at 1300 km with decay pipe 2m radius 180 m length 


0° 


52.5 


19.4 


0.374 


0.074 


3.87 


1.05 


0.223 



9. DETECTOR REQUIREMENTS 

The detector requirements for a detector in a beam towards DUSEL and a detector in the NuMI 
off-axis beam are quite different. Although the physics goal of measuring 613, mass hierarchy, and, 
above all, CP violation is the same, the obstacles to obtain sufficient sensitivity to this physics are 
very different for the two techniques. We will describe the understanding reached in the process of 
this study. 

Both techniques are attempting to obtain sensitivity to CP violation in the neutrino sector by 
collecting sufficient numbers of Vg appearance events. By obtaining appearance events at 

difference oscillation phases and energy, matter effects and CP effects can be disentangled to mea- 
sure oscillation parameters without correlations or ambiguities. Regardless of the technique the 
most important experimental parameters are the numbers of events at or near the oscillation peaks 
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TABLE VII: Signal and background anti-neutrino interaction rates at 1300 Km (Fermilab-HOMESTAKE) 
using different WBLE beam energies and off- axis angles. The rates integrated over the neutrino energy 
range of - 20 GeV. Rates are given per MW.lO^s.kT. For v^^ — > Vg oscillations a value of sin^20i3 = 0.04 
and Arn^y = 2.5 x 10^^ eV^ is used. No detector model is used. 



Degrees off-axis 


v^CC 


CC osc 


Ve CC beam 


Vg QE beam 


NC-l7r° 


v^, ^ Vg CC 


Ve QE 


WBLE 120 GeV at 1300 km with decay pipe 2m radius 380 m length 


0° 


75.0 


37.7 


0.570 


0.106 


7.79 


0.669 


0.160 


0.5° 


33.5 


13.0 


0.356 


0.077 


4.90 


0.332 


0.103 


1.0° 


12.0 


6.47 


0.185 


0.056 


2.64 


0.122 


0.056 


2.5° 


1.41 


0.694 


0.037 


0.013 


0.499 


0.033 


0.022 


WBLE 60 GeV at 1300 km with decay pipe 2m radius 380 m length 


0° 


50.5 


21.3 


0.373 


0.088 


6.05 


0.507 


0.137 


0.5° 


25.4 


8.52 


0.248 


0.066 


4.23 


0.272 


0.094 


1.0° 


10.3 


5.38 


0.144 


0.045 


2.52 


0.116 


0.058 


2.5° 


1.36 


0.667 


0.031 


0.013 


0.518 


0.035 


0.024 


WBLE 40 GeV at 1300 km with decay pipe 2m radius 380 m length 


0° 


33.8 


12.5 


0.270 


0.069 


4.70 


0.366 


0.110 


WBLE 28 GeV at 1300 km with decay pipe 2m radius 180 m length 


0° 


14.6 


4.94 


0.076 


0.026 


2.64 


0.172 


0.065 



versus the numbers of irreducible and reducible backgrounds. The numbers of events in either 
technique are roughly proportional to the exposure defined as the beam power in MW (at some 
chosen proton energy) times the total detector efficient fiducial size in kTon times the running time 
in units of 10^ sec. In the following, to set the rough scale for detectors, we will assume that a few 
hundred events after accounting for detector efficiency are needed at sin^20i3 = 0.1 per 

year. As pointed out in Section [5| accelerator power of ~ 1 MW can be obtained and handled with 
current technology; this sets the scale for the detector size, efficiency, and running times. 

9.1. Off-axis 

In the off- axis technique, we have considered two large detectors at two different locations. On 
the NuMI beam-line, the places considered for the placement of these detectors are: 1) baseline 
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length of 8 10 km and off-axis distance of 12 km, 2) baseline length of 8 10 km and off-axis distance 
of 40 km. At a length of 810 km (which is close to the maximum possible on the NuMI baseline), 
the first and second oscillation maxima for the physics under consideration are at neutrino energy 
of 1.64 GeV and 0.54 GeV, respectively, for 5m^2 = O.OOlSeV^. The off- axis distances were 
chosen to obtain a narrow band neutrino beam at or near these oscillation maxima. These spectra 
and the event rates can be seen in [fT2|. 

Shorter baseline lengths for NuMI off- axis detectors have been considered in the literature [|34l . 
We have commented on this approach as part of the answers to questions in Appendix A. We will 
not consider this approach here because of the practical difficulties noted. 

The main detector requirements for off-axis detectors are: 

• Size: To approach the exposure criteria of few hundred events per year for sin^ 20i3 = 0. 1 the 
total efficient fiducial mass of the detectors at the first and second oscillation maxima needs 
to be ~ 100 kT. This could be deployed with 50 kT at the first location (12 km off-axis) and 
50kT at the second location (40km off-axis) or all of the mass in one location. 

• Cosmic ray rejection: NuMI based off-axis detectors will likely be on the surface or have 
a small amount of overburden. Surface or near-surface capability is essential for the NuMI 
based off-axis program because of the geographic nature of the area. As pointed out in 



Section 12 a surface detector needs to a) have sufficient data acquisition bandwidth to collect 
all events near the beam spill time, b) eliminate cosmic ray tracks so that the beam events can 
remain pure, c) tag events due to cosmic rays so that no cosmic ray induced events mimic an 
in-time beam event. These requirements force the surface detector to be a highly segmented 
detector with active cosmic ray veto shielding. 

Background rejection: There are two contributions to the background from the neutrino 
beam: neutral current events and contamination of electron neutrino events. The narrow 
band nature of the neutrino beam is important for rejection of both of these backgrounds. 
The neutral current events which tend to have a falling energy distribution can come from 
both the main peak of the neutrino spectrum and the tails. In the case of the second location, 
40 km off-axis, the large kaon peak will contribute background. The Ve contamination has 
a broad distribution for both off-axis locations [|T2ll . To use the narrow band nature of the 
beam effectively to suppress backgrounds, the detector must have the capability to measure 
neutrino energy (total charged current event energy) with good resolution, which is approxi- 
mately the same as the width of the narrow band beam. It should also be able to reject or 
photon induced showers. 
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9.2. Detectors at DUSEL 



The two sites for DUSEL that made a presentation to this study are 1290 (Homestake) and 
1495 (Henderson) km from FNAL. The study has considered distances as far as ~2500 km and 
concluded that the physics capability, with some exceptions, is roughly the same for same sized 
detector. The first and second oscillation maxima for 1290 km are at 2.6 GeV, and 0.87 GeV; 
for 1495 km, they are at 3.0 GeV and 1.0 GeV, for Am^2 = 0.0025eV^. A new neutrino beam 
at O'' or at small off-axis angles has been simulated liT4l to show that a spectrum could be made 
to cover these energies; the critical parameter in the flux at low energies will be the decay tunnel 
diameter which must be kept to be ~ 3 — 4m, which is a factor of 1 .5-2 larger than the NuMI decay 
tunnel. The beam-line could be operated at any energy between 30 to 120 GeV proton energy. For 
higher proton energies work is in progress to remove high energy neutrinos (> AGeV) that produce 
background. The beam-line could also be operated at a slight off-axis angle if the background can 
be lowered by modest amount while operating at the highest power level possible at 120 GeV. For 
the purposes of setting broad detector requirements we will assume that the spectrum is similar to 
Figures |3] or [6} 

Detectors at DUSEL (at either Homestake or Henderson) could be placed either on the surface 
or at a deep site. If placed on the surface the detector considerations would be approximately the 
same as those for off-axis detectors because the primary design issue would be rejection of cosmic 
ray background. The availability of deep sites at the appropriate baseline distance for a very large 
detector are the main reason for locating the detector at DUSEL. Both Henderson and Homestake 
are planning on large detector caverns at a depth of ~ 5000 ft. We will enumerate the detector 
requirements assuming this depth. 

• Size: To approach the exposure criteria of a few hundred Vg appearance events per 
year at sin 20i3 = 0.1, the efficient fiducial mass of the detector needs to be ~ IQQkT . In 
the case of DUSEL all of this mass can be in the same place exposed to a beam that contains 
both oscillation maxima. 

• Cavern: Because of the size required for the detectors, a stable large cavity (or cavities) 
that can house ~ \QOkT of efficient fiducial mass will be needed. For a water Cherenkov 
detector, which is well suited for deep operation, the efficiency is expected to be ~ 20% 
indicating a real detector size of several hundred kTon. From preliminary studies it appears 
that both Henderson and Homestake satisfy this criteria. 

• Cosmic ray rejection: Since the cosmic ray rate at the deep sites proposed for DUSEL detec- 
tors is very low, it will not be a major factor in detector design. A cosmic ray veto for such 
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a detector might be needed for physics other than accelerator neutrino physics; for example, 
detection of solar neutrinos. But it is not required for the physics discussed here. 

• Surface location for a detector: For a liquid argon TPC, the efficiency and background re- 
jection could be high and therefore the detector could be ~ lOO^r. However, for an un- 
derground liquid argon TPC the requirements on the cavern will be dominated by safety 
concerns regarding storage of such a large amount of cryogenic liquid in a deep laboratory. 
If the liquid argon detector is placed on the surface, the requirements are approximately the 
same as for the NuMI based off-axis detectors. The dominant requirement will be rejection 
of cosmic ray background. 

• Background rejection: There are two main contributions to the in-time background from 
the beam: neutral current events, and electron neutrino contamination in the beam. It is 
expected that the majority of the NC background at low energies will be from single Tt^ 
events that will have to be rejected. In the case of using a wide band beam, there are two 
tools for signal extraction. Pattern recognition with good capability will be needed to reduce 
neutral currents, especially single Tt^ events. The oscillation pattern in the energy spectrum 
will also be used to extract the signal. The first oscillation node, in particular, will form a 
peak above 2 GeV with a well known shape. To allow such a signal extraction, the detector 
must have good energy resolution for neutrino energy. From the work reported here ~ 10% 
energy resolution above 0.5 GeV including Fermi motion effects will be needed. For a water 
Cherenkov detector there is new work on pattern recognition to reduce the NC backgrounds 
and obtain the needed energy resolution. For a liquid argon detector, it has been shown that 
the NC background can be suppressed to very low levels for low multiplicity events (such as 
quasi-elastics) while maintaining good resolution. 

10. STATUS OF DETECTOR SIMULATIONS 
10.1. Water Cherenkov Detector 

As part of this work, we have studied the background rejection and neutrino energy resolution 
(from charged current events) of a large water Cherenkov detector instrumented in the same manner 
as Super-Kamiokande. Although considerable further work is needed the capabilities appear to be 
sufficient for the neutrino oscillation program under consideration. The total mass and exposure 
needed to achieve good sensitivity to CP violation in neutrino oscillations was also determined. 

The technique of water Cherenkov detectors with non-focusing optics is well understood. In 
particular, the light yield and the fraction of scattered light can be modeled accurately. Software 
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techniques exist that use the pattern of light and the time sequence of photons to reconstruct vertices 
and trajectories of charged particles. The vertex resolution depends on the timing accuracy of the 
PMTs. The energy resolution and the energy threshold depends on the total amount of detected 
hght. Both of these have been extensively discussed in technical articles and Ph.D. theses [|35l[36l . 
Considering the substantial existing knowledge and information about this technology, we decided 
to focus only on the additional new requirements imposed by the accelerator neutrino physics 
under consideration. 

For the program considered here an essential problem is to separate electron shower events 
from neutral current events, especially events containing a single in the final state. The goal is 
to search for Ve charged current induced showering events in the 0.5 to 4 GeV range. For example, 
single ^r'^ particles with energies of 1, 2, 3, and 4 GeV decay to two photons with a minimum 
(which is also the most probable) opening angle of 16, 8, 5, and 4 degrees, respectively. The 
probability of a decay with an opening angle of more than 20° for 1, 2, 3, and 4 GeV ;r^'s is 
40%, 8.2%, 3.6%, and 2.0%, respectively. In a water Cherenkov detector the position where the 
7r° photons convert cannot be measured with sufficient precision from the pattern of Cherenkov 
light, which tends to be two overlapping showering rings. At low energies the opening angle is 
sufficiently large compared to the Cherenkov angle (~ 42°) that single ;r^'s can be separated quite 
effectively. At energies greater than 2 GeV, however, the small angular separation between the two 
photons makes such separation difficult. It is well known that resonant single pion production in 
neutrino reactions has a rapidly falling cross section as a function of momentum transfer, q^, up 
to the kinematically allowed value ll37l . This characteristic alone suppresses the background by 
more than 2 orders of magnitude for n'^ (or shower) energies above 2 GeV. Therefore a modest 
background suppression (by a factor of ~ 15 below 2 GeV and ~ 2 above 2 GeV) should make the 

background manageable over the entire spectrum. 

As part of this study such background suppression has been demonstrated using complete simu- 
lation and reconstruction using the Super- Kamiokande detector as the benchmark [|38l[39ll . Similar 
suppression has also been obtained independently by another group [|40l |4T|| . In both studies the 
rejection of backgrounds was enhanced beyond the currently well known capabilities of a Super- 
Kamiokande like detector by using a combined likelihood method. In this method a number of 
event observables (a complete list can be obtained from the talks in [|38l and [|40l ) with low back- 
ground discriminating power were combined in a single likelihood cut. The work in [|38l chose to 
cut on this likelihood as a function of reconstructed energy so that the efficiency of this additional 
cut for charged current electron neutrino events was ~40%. The additional rejection for neutral 
current events ranges from a factor of 30 at 300 MeV to a factor of 4 at 3 GeV (page 34 in [|38ll ). 



Table VIII shows the rejection power achieved by this method as a function of energy. The table is 
divided in two parts: before the event energy can be reconstructed the rejection can only be given 
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in terms of true quantities such as "true energy". After the event is reconstructed the rejection is 
given in terms of "reconstructed energy". It should be remarked that the reconstructed energy for 
an NC event is considerably lower than the true neutrino energy. The total integrated efficiency for 
signal in this calculation using the 28 GeV spectrum is 37% from the traditional cuts multiplied by 
40% efficiency of the likelihood cut. As explained in [|38ll the likelihood cut could be adjusted to 
have higher efficiency at a cost of higher background. The total integrated rejection of neutral cur- 
rent background is ~ 13 for the traditional cuts multiplied by 24.0 using the subsequent likelihood 
cutsll42l. 

The work reported in POl compares her results to [|38l (page 31-31 (4U\). In the work reported in 
[|40ll it was chosen to retain high efficiency (above 70%) to electron neutrino charged current events; 
she obtained rejection factors of 10 at 300 MeV declining to 2 at 3 GeV. The two calculations are 
in good agreement if compared at the same efficiency considering that the simulated event sample 
and methods of discrimination were quite different. 

The other important component of this study is the neutrino energy resolution for charged cur- 
rent electron neutrino events. The selection procedure described in the previous paragraphs at- 
tempts to select clean events with a single lepton in the final state. To measure the neutrino energy 
we assume that this selected event is a quasi-elastic scattering event. We then calculate the neutrino 
energy using the following formula: 

'^i^p E lepton "t" Plepton ^lepton) 

to reconstruct the neutrino energy (fi^ec) from the measured electron energy (E[gpfon) and electron 
angle (Oiepton) with respect to the neutrino direction. The energy resolution using this method has 
four components. The energy resolution of the electron has been demonstrated to be 4% at 500 
MeV improving to 2 % above 2 GeV (page 84 in [[36ll '). The angle of the electron with respect 
to the beam must be measured to calculate the energy of the incoming neutrino. The angular 
resolution ranges from 3 deg at low energies to 1.5 deg at high energies (page 81 in [|36ll '). The 
third component to the neutrino energy resolution is the Fermi motion of the struck nucleon inside 
the oxygen nucleus. This is often modeled using either data from electron scattering or using a 
simple Fermi gas model. It adds a contribution of ~ 100-200 MeV to the resolution. Finally, the 
selected events have a contamination of non-quasielastic events in which the extra particles (such as 
charged mesons or photons) in the final state are either invisible because they are below Cherenkov 
threshold or are missed because of poor reconstruction. The final energy resolution including all 
these effects has been calculated (page 17 in P3ll ) to be about ~ 10% at 1 GeV with significant 
non-gaussian and asymmetric tails. These tails are due to the nuclear effects and non-quasielastic 
contributions. The resolution improves at higher energies. The effect of the resolution is that the 
oscillation pattern remains visible although somewhat degrades. The resolution will have to be 
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Cut 


Energy Bin (GeV) 




0-0.5 


0.5-1.0 


1.0-1.5 


1.5-2.0 


2.0-3.0 


>3 


SK cuts/ E,rue 














Ve signal 


74% 


74% 


62% 


44% 


36% 


27% 


CCv^ bkg. 


0.17% 


0.44% 


0.75% 


0.76% 


0.90% 


0.45% 


NC bckg. 


1.7% 


4.43% 


5.3% 


7.0% 


7.7% 


8.6% 


Beam Ve bkg. 


86% 


75% 


57% 


46% 


36% 


23% 


Likelihood cuts/^rec 














Ve signal 


40% 


40% 


40% 


39% 


40% 


40% 


CCv^ bkg. 


6.8% 


13.6% 


6.3% 


8.0% 


6.5% 


2.2% 


NC bkg. 


0.72% 


4.5% 


6.3% 


3.9% 


8.3% 


7.0% 


Beam Vg bkg. 


37% 


41% 


40% 


37% 


39% 


34% 



TABLE VIII: Simulation and analysis results on the fraction of events kept after the traditional cuts (top 
part of the table) and the additional efficiency after the newly developed likelihood cuts (bottom part of 
the table). The events are divided in 4 parts: signal from Vg charged current events (of which a small 
part are quasi-elastics), charged current events, neutral current (NC) events, and background due to Ve 
contamination in the beam. There is no entry for background from charged current because the beam 
spectrum is dominantly below T production threshold (~ 3.5GeV) and this background is estimated to be 
low. The efficiency for the signal and Ve background should be the same except for the statistical fluctuations 
in the Monte Carlo due to small statistics of the Ve background. We have retained the numbers in the table 
to demonstrate consistency. 

modeled well to extract the oscillation signal and the oscillation parameters from the far detector 
data; this is true regardless of the detector type. 

The sensitivity calculations described later in this report were performed using the GLoBES 
framework [|44l . For this calculation the detector response was parameterized and adjusted to cor- 
respond to the full simulation described above. There are some differences that should be kept 
in mind to allow comparison between calculations. The first difference is that in the sensitivity 
calculation no events below 0.5 GeV are used. The work in J38l and [|40l includes events to lower 
energies. The second difference is in the energy resolution. The energy resolution obtained after 
complete simulation and reconstruction is shown in Figure [8] bottom plot. The sensitivity calcula- 
tion has a parameterized resolution function that includes effects of Fermi motion, resolution on 
the lepton energy and angle, and non-quasielastic contamination. The two resolution functions 
are shown in Figure [8] The parameterization has somewhat worse resolution in the core than the 
full simulation, but less tail than the full simulation. The input to the calculation is firmly based 
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on full simulation. But the parameterized background and resolution allows for fast calculation of 
signal and background rates for different oscillation scenarios; we are also able to change the beam 
spectrum while keeping the detector performance fixed using this tool. 

We note that comparison of signal rates between various authors needs to be done with some 
care because they may have used different conventions for the sign of the CP violating phase 
dep. The calculation here uses the standard convention from PDG (with e in the t/gs matrix 
element), but some of the simulations reported in the course of this work have the convention with 
the opposite sign on 5. 

We have explored the feasibility of using a water Cherenkov detector for this science with 
promising results, but considerable further work is needed to optimize the detector and gain com- 
plete understanding of its limitations. One question that needs exploration is what are the intrinsic 
hmitations of the water Cherenkov technique in terms of pattern recognition. The literature on the 
subject is broad and general conclusions can be drawn [|45l l46l l47l |48l |49l. The question can be 
quantified in terms of the vertex and angular resolution for single tracks and the ability to separate 
two tracks that are close in angle and have a common vertex. In the references cited above it has 
been shown that the photon detection resolution (time and position) in a Super-Kamikande style 
water Cherenkov detector (proximity focused) does not approach the multiple scattering contri- 
bution to that resolution. Therefore, the capabilities of the detector could be improved by modest 
improvements to the timing and granularity of the PMTs or by addition of ring imaging techniques. 
In other words, the current capabilities are dominated by the characteristics and geometry of the 
photo-multiplier array, and there is room for improvement. 

Photomultiplier coverage and optimization: The optimization of PMT coverage and granu- 
larity has not been addressed in our study. This is an appropriate goal for the proponents of the 
water Cherenkov technique when they write the full proposal. Some hints of the effects of high 
granularity can, however, be found in the work of Yanagisawa and Dufour. For example, in [l38ll it 
is shown that the efficiency for detecting a tjP particle increases by 20% as the event moves away 
from the wall of the detector; this indicates that a larger detector or a detector with more granular- 
ity with the same PMT coverage will have better background suppression. In [|50l . the likelihood 
based background suppression is shown to have weak dependence on the PMT coverage (either 
20% or 40%). These preliminary results indicate that as long as the collected numbers of photo- 
electrons is reasonably large, the granularity of the PMT readout will have more impact on pattern 
recognition. 

Optimizing beam spectra: In this study we have not made extensive attempts at optimizing 
the beam spectra versus the detector performance for backgrounds. In [fT4ll . a cross section model 
was used to calculate signal and background shapes for neutrino spectra produced with different 
energy proton beams. It is clear that the background will increase for proton energies above 60 
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FIG. 8: Top is the parameterized resolution function used in the calculation of sensitivity. Bottom is the 
energy resolution on selected electron neutrino events after all cuts as described in |[39ll43l . The tail where 
reconstructed energy is lower than true energy is due to non-quasielastic events that are selected. These 
events have missing particles and therefore have missing energy. The bottom plot was made from simulated 
events including the effects of oscillations over 1480 km, and therefore there is a depletion of events in the 
plot around the oscillation minimum. 
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GeV because of the high energy tails. However, for proton energies below 60 GeV, the signal 
versus background performance is approximately constant: there is reduction of background for 
lower energies, but there is significant statistical gain (for constant power) at higher energies. The 
intensity at low neutrino energies as well as the elimination of high energy tails will continue to 
need examination. 

There are three ways to optimize the beam spectrum to reduce backgrounds from high energy 
neutrinos: 1) The optics of the target horn assembly need to be optimized to increase the flux 
below 3 GeV and reduce it at higher energy, 2) The beam towards DUSEL could certainly be made 
slightly off-axis. Reference [[T4l shows that a 0.5" off-axis angle lowers the event rate > 5 GeV by 
about a factor of 3. The 0.5° off-axis beam angle can be accommodated in a large 4 meter diameter 
tunnel, but the implications for shielding, beam dump, as well as future flexibility of the program 
must be carefully considered. 3) The third option for reducing neutrino flux at higher energies is 
to introduce a second target (or beam plug) between the two focusing horns. This has been studied 
for NuMI and found to be effective at reducing tail events by as much as 70% [ISTl . Such an option 
needs careful engineering design because it will affect the radiation environment in the target area. 

Near Detector : Lastly, there has been very little study of near detector issues for a beam to- 
wards DUSEL. There are 3 issues that need to be discussed. 1) The availability of space: the 



beam design discussed in Section 7.2 leaves 300 m of earth shielding from the end of the decay 
tunnel to a potential near detector site located within FNAL boundaries. The depth at this site 
would be 192 m for a beam to Henderson and 176 m for Homestake. The NuMI near detector 
is at a depth of 105 meters. The feasibility and cost of a near detector cavern of about 30 meter 
width/length will need to be examined. 2) The main requirement for the near detector: The most 
important function of the near detector is the measurement of the neutrino spectrum and back- 
grounds before oscillations. As explained later for the sensitivity calculation we have assumed 
that the background will be known to about 10%. This includes the effects of beam and detec- 
tor modeling as well as nuclear effects which might be different between the near and the far 
detectors if they are composed of different materials. A harsher requirement on the near versus 
far energy scale of < 1% might come from the need to measure Am^j with high precision. 3) 
The detector granularity, mass, and data acquisition: the event rate at the near site (few events per 
beam pulse for a ~ 100 ton detector) will be much higher compared to the far detector. Reconstruc- 
tion of these events will likely require a fine grain detector (perhaps a modest sized hquid argon 
TPC) with electronics that can separate events within the 10 lis pulse. How the design impacts the 
requirements outlined above will need to be examined. There is now considerable experience from 
the NuMLMINOS on how to use the near detector to perform the appropriate extrapolations from 
near to far site. From that experience [|7l the above requirements do not appear to be particularly 
difficult, but the issue should not be treated cavalierly. 
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10.2. Liquid Argon Time Projection Chamber 



Liquid Argon TPC detectors, with fine grained tracking and total absorption calorimetry capa- 
bility, suggest great promise for sensitivity to long baseline oscillation physics. Hand scanning 
studies indicate efficiency for charged current quasi-elastic electron neutrino interactions (Vg QE) 
greater than 80% and background rejection of neutral current ^r*^ events by a factor of 70 [52] . 
Studies from European groups are consistent with these results [|53l . 

As part of this study, tools have been developed to simulate and reconstruct events in the Ey = 
0.5-5 GeV energy region. Studies using these new tools confirm the efficiencies and background 
rejection from the hand scanning work. Sensitivity calculations folding in these efficiencies, back- 
ground rejection factors, and resolutions indicate LArTPCs are ~3 times more sensitive than an 



equal mass of Water Cherenkov detector (See Section 13). 

It is primarily the imaging capability that enables LArTPCs to distinguish different event classes 
from each other. Specifically, while conventional detectors can typically identify only the outgo- 
ing lepton in QE interactions, LArTPCs can tag both the outgoing lepton and the recoil proton. 
Furthermore, a LArTPC can easily and unambiguously identify the interaction point of energetic 
gamma-rays, for example from ;r° decay, if the separation from the primary vertex is larger than 2 
cm|l54l. 

As part of this study, a GEANT3 simulation of a Liquid Argon TPC was studied and developed 
to best quantify the detector performance [|55l . The Monte Carlo used the NUANCE event genera- 
tor as input and simulated events in a 7 x 10 x lOm^ box, roughly equivalent to 1 kTon. Events are 
digitized using standard GEANT libraries, and Monte Carlo truth studies performed on this out- 
put. Given the imaging capability of a LArTPC, this is an acceptable approximation of an actual 
event. The criteria to tag a Vg QE interaction are first to see an electron shower as distinct from 
a muon track. This is assumed to be 100% efficient. The second criteria is to see a recoil proton 
coming from the same vertex as the electron. The well established low energy threshold for this is 
a proton with kinetic energy > 40 MeV lf56l . Imposing this requirement, the efficiency for QE 
events is >90%. As these are first pass studies, we default to the more conservative 80% efficiency 
determination from the hand scanning. Neutral current backgrounds, with subsequent n*^ 77, 
arise from both v^n — > Vfj,nn^ and v^p — » v^pn^ interactions. The first, V^n v^nn^, is rejected 
because of the lack of any recoil proton. The second, v^ip v^ipK^, is tagged by observation of a 
2cm or larger gap between the vertex of the recoil proton and at least one of the gammas from the 
decaying which converts into an e+e^ shower. Combining these requirements, only 0.5% of 
the NC ^r*^ backgrounds are not rejected. Further rejection factors are expected by looking at the 
energy deposited in the first few cm of tracks initiating electron showers versus gamma showers. 
The overlapping e+e^ from the gamma shower deposits twice the energy at the beginning of the 



43 



track as the single electron. 

These studies have been performed using the WBLE flux generated with 40 GeV protons used 
in this study. Further study is needed to understand reconstruction for the other beam options and 
the NuMI option. Nevertheless, these results are relevant across a broad range of energies. In 
particular, for NC rejection, the separation between the primary vertex and the closest gamma 
conversion point is roughly independent of the incoming neutrino energy [|54l . High multiplicity 
events in the deep inelastic scattering (DIS) region may be very challenging to reconstruct. Effi- 
ciency and background rejection for DIS events for the different flux configurations is also needed. 

Advances in automated reconstruction were also pursued as part of this study. The Hough 
transform based fit algorithm was designed to reconstruct linear tracks from a quasi-elastic event 
through a parameterization by angle. It efficiently identifies both primary and secondary vertices 
and reconstructs tracks with resolution of ~ 2° (RMS) [|54l 1571 . This fitter suffices for events 
with linear tracks and low multiplicity such as quasi-elastics and resonance events. A study of the 
capability to automatically identify and reconstruct electromagnetic showers is in the early stages. 

In the future, this simulation and reconstruction package can be used to study energy resolution 
for different classes of events. 

For the results in this report, we use the energy resolution from previous work. For QE events, 
a 5% energy resolution was assumed. This is valid down to ~1 GeV, below which few events 
contribute to the oscillation signals. For non-QE events, a 20% neutrino energy resolution was 
assumed. This is likely too conservative in the resonance region where low multiplicity events can 
still be well measured by LArTPCs, but likely too optimistic for DIS events above 2-3 GeV. Un- 
derstanding these resolutions as a function of energy is part of the ongoing program of simulation 
and reconstruction studies. 

11. STATUS OF DETECTOR DESIGN AND TECHNOLOGY 
11.1. Water Cherenkov conceptual Design 

The water Cherenkov detectors discussed in this study were largely conventional based on the 
well known technology developed and perfected over the last three decades. The main difference 
is the factor of ~ 10 increase in fiducial mass compared to the largest existing detector (Super- 
Kamiokande). This large increase can be accomplished either by increasing the size of the detector 
or by building several detectors (or both). The second important parameter for this detector is the 
number and size of the photo-multiplier tubes (PMTs). 

Two conceptual designs were reported for this study. They were specifically for the two possible 
DUSEL locations of Homestake or Henderson, but the authors have acknowledged that their ideas 
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could be adapted to either site with appropriate considerations for site dependent cost factors. 

The design reported for the Henderson site (UNO [|58l ) has a single cavity of dimension 60 m 
wide, 60 m high, and 180 m long. The 180 meter length is divided in 3 sections. Each section is 
a separate optical volume with photo-cathode coverage of 10% for the end sections and 40% for 
the central section using the 20 inch diameter PMT developed by Hamamatsu. Each section has 
fiducial volume (depending on specific physics cuts) of about 150 kTon. The depth of the detector 
in Henderson will be approximately 5000 ft. 

The design reported for Homestake houses the detector at 4850 ft depth in 3 separate large 
caverns or modules [[59ll . The size of the caverns will be cylindrical with diameter/height of ~ 
53 meters. The location is at the 4850 ft level of Homestake which is proposed as the Early 
Implementation Plan for the Homestake lab. The collaboration proposes that the same level be used 
to accommodate several more cavities to take the total detector mass to megaton over a long period, 
but the baseline detector is 3 modules. The exact dimensions of the cavities will be determined by 
the need to maintain fiducial mass of 100 kTon for accelerator neutrino events. Each detector will 
be instrumented by 10 to 13 inch diameter PMTs with photo-cathode coverage of 25%. At this 
stage of simulation and understanding of PMT performance, the Homestake proponents consider 
the choice of smaller but larger numbers of tubes for granularity adequate for reconstructing the 
accelerator neutrino events. The concept for the Homestake detector including the physics and 
a rough estimate for the cost was presented to a program committee for the Homestake interim 
laboratory. The review can be obtained at [|60l . 

The main concerns for both designs is the cost and time required to build stable and safe 
cavern(s)[|61j and the manufacturing of the necessary number of photo-multiplier tubes. 

For the single cavity (UNO) concept an estimate based on the cost of Super- Kamiokande has 
been made for the cavity excavation of $168M; the engineering and stability of the cavity needs 
detailed examination. The total cost including 56000 large 20 inch PMTs and 15000 smaller 8 
inch tubes for outer veto volume was estimated to be $437M. The total construction time will be 
approximately 10 yrs dominated by the PMT manufacturing time[58|. 

For the multi-cavity Homestake design, the proponents have performed an initial engineering 
design for the cavity construction and a stability study [|62l . The cost for constructing 3 cavities is 
estimated to be approximately $70M which includes contingency factors. The time scale for con- 
structing the first cavern is 4 yrs and each additional cavern is readied 6 months after the completion 
of the previous one. The total cost including approximately 50000 PMTs for each detector module 
is $309M. The impact on this cost if the size of the module is increased for additional fiducial 
volume is explained in Appendix A. Based on the Super-Kamiokande experience, the installation 
time for the PMTs will be about 1 yr for each module. 

The largest unknown at present for both designs is the schedule for manufacturing the large 
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numbers of PMTs. For the 20 inch PMT option, there appears to be only one vendor at present 
with a labor intensive manual process. For the smaller PMTs there could be multiple vendors with 
more automated manufacturing processes. We will comment on this issue again at the end of the 
report. 

There are other technical concerns for such a large water Cherenkov detector: the handling, 
temperature and purification of such a large amount of water, the engineering for mounting the 
PMTs and the cabling of the large number of channels, maintenance of the PMTs and associated 
electronics, and the radiation environment in the deep site which can affect the data rate and the 
energy threshold of the detector. There is no detailed engineering design for these items, however 
these issues have been examined by previous generations of these types of detectors. Based on that 
previous experience both detector designs have included approximate costs in their estimates. 

11.2. Liquid Argon TPC Conceptual Design 

While LArTPCs show great promise with excellent efficiencies and background rejection for a 
variety of physics goals, they have not yet been demonstrated on scales larger than few hundred 
tons in size. An active R&D program culminating in the T600 program [|63l has illustrated the 
capabihties of the detector, however, further R&D is necessary to consider massive detectors, on 
the scale of tens of ktons. 

There are several different design ideas for massive detectors including a modularized detec- 
tor [[64l . a single detector but with modularized drift regions [52], and a single open volume, very 
long drift detector combining charge and light collection [53 1. The technical issues described here 
are relevant primarily for single massive detectors with modularized drift regions, the design stud- 
ied by the contributors to this study. For these, there are no major obstacles to scaling to detectors 
on the scale of 50-100 kTon, however, there is an R&D path that must be realized in order to con- 
sider massive detector construction, operation, and data analysis. Details of this path and major 
R&D goals can be found in f,65.l . The major challenges for scaling to a large detector include: 

• Argon purity 

• Signal to noise in a massive TPC 

• Understanding Cost and Schedule 

Progress and path for each of these is described below. 

For ionization electrons to drift 3 m in a LArTPC, 10ms electron lifetime must be achieved and 
maintained. Studies from the T600 run suggest this is possible, but, for a massive detector, mod- 
ifications must be made to the purification system, and the ability to reach purity levels necessary 
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in an industrial environment must be demonstrated. Over the past year, Fermilab has embarked 
on purity testing towards this goal. They have developed a new non-proprietary Trigon filter (un- 
hke those used in the T600) that can be regenerated in-line. With this filter system, Fermilab has 
achieved 12 ms lifetimes in a small test vessel. Over the next year, purity studies will continue with 
a materials test stand [f66ll at Fermilab where argon will be re-purified after being exposed to con- 
taminants expected in a massive LArTPC. An additional challenge to purity is a consequence of 
the inability to achieve vacuum before the initial argon fill in a massive detector. An idea to purge 
the vessel with clean argon gas prior to liquid fill is being tested at Fermilab now with studies 
continuing in the upcoming year [67|. 

A very massive detector will have signal wires as long as tens of meters. Long wires present 
challenges related to wire breakage, wire assembly and stringing, and electronics noise. Existing 
R&D work at Fermilab focuses on assembly techniques and noise pickup using a long wire test 
stand ll68l . Work on electronics design to maximize signal to noise specifically by employing cold 
electronics, is underway at Fermilab and Michigan State University. A new idea for internal wire 
configuration, a cellular design, avoids many of the stringing and assembly problems of long wires 
by stringing wires onto pre-assembled ladders before installation [[69ll . 

There are two cost drivers for a liquid argon TPC which have some certainty at this point. The 
first of these was given by the LArTPC group in its September 2005 report to NuSAG [|52l . There, 
the cost of liquid argon alone (without a purification system) is reported as about $1 M per kton. 
Subsequent to that report, a simple scaling relationship has been developed based on information 
from two vendors for tanks appropriate for containing liquid argon (but without modifications 
required to put a TPC inside it). This relationship, which is expected to be valid between 5 kTon 
and 50 kTon, is $2.72M -i- 0.306 $M/kTon. Thus taking a 50 kTon detector as an example, these two 
cost drivers (the liquid argon plus a containment tank) would cost about $50M + $18M = $68A/. 

There are many other costs, both technically driven and project driven, but the design of the TPC 
itself needs to be specified in more detail before such a complete costing exercise can converge. For 
example, the recently developed cellular design for the TPC significantly changes the requirements 
on the containment tank compared to the design in the September 2005 LArTPC report to NuSAG. 
Since this design allows for fabrication of the TPC wire planes at the same time as the containment 
tank is being constructed, the schedule for construction of the detector is shorter. If electronics are 
used in the liquid argon, the cellular design will change and the requirements on pattern recognition 
will become easier. Finally, the idea of using several smaller tanks to achieve a large mass will 
impact the cost of the purification system as well as the cost of the containment tank(s). These are 
some of the design choices for the TPC that need to be made before a cost estimate of the technical 
components, other than the liquid argon and the cost for a single containment tank, can be made. 

In addition to the major challenges for scaling to large detectors as described above, issues 
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relating to detector siting have been studied. Water Cherenkov detectors must be located deep 
underground due to cosmic ray backgrounds. By contrast, liquid argon detectors could be located 
on or near the surface. As part of this study, cosmic rates in a massive LArTPC detector were 
calculated and their impact on the physics program was considered [TTOII and is discussed in more 



detail in Section [T2j If massive LArTPCs are sited with some overburden, such as at the 300 ft 
drive-in site at Homestake, cavern construction for these detectors must still be understood. As 
part of this study, cavern designs modeled after liquefied natural gas vessels built within ships 
hulls were considered [TtTII . This design is promising and studies on this are ongoing. 

The R&D path towards a massive detector includes small scale tests and studies as described 
above. Construction of a significantly larger prototype, ~lkTon, is necessary before embarking on 
the massive detector project. The details of this R&D path at Fermilab will be addressed within 
the next year. 



12. OVERBURDEN AND SHIELDING 

In this section we briefly discuss the overburden issue in the context of accelerator neutrinos. 



For non- accelerator physics the issue is discussed in Section 14 

In summary, the background rates in a large detector due to cosmic rays have been calculated for 
both surface and underground locations. A preliminary evaluation of the consequences for both 
data acquisition and background to accelerator neutrino events suggests: 1) It is not possible to 
operate a water Cherenkov detector of size > 50kT on the surface. 2) A fine grained tracking 
detector such as a liquid argon TPC could be operated on the surface to take data within the 
short (~ 10/i5 at FNAL Main Injector, ~ l.SjlS for BNL AGS) accelerator spill[|70l|, however 
background rejection of ~ 10^ (~ 10^ — 10^) will be needed against cosmic muons (photons) by 
either active veto or pattern recognition to reduce the background rate to acceptable levels; this 
rejection is in addition to the rejection obtained by the timing requirement. We provide a few more 
details of the calculations below. 

A cylindrical tank of size 50 m height/diameter (approximately lOOkT of water) will have a 
rate of cosmic muons (with momentum > O.SGeV /c) 250 kHz from the top and 250 kHz from 
the sides. For a 10 jls beam spill this corresponds to 5 muon tracks in the detector. For a single 
volume water Cherenkov detector in which the photo-multipliers are mounted on the walls looking 
inwards, each muon on the average will produce a hit in more than 50% of the PMTs. Therefore, 
each cosmic ray should be assumed to deaden the entire detector for a period of time which is 
dependent on the dwell time of the muon track and the light inside the detector, the pulse shapes 
from PMTs, and the data acquisition electronics. All the effects are of order IjAs and therefore 
make the detector unworkable at the surface. For example, for a detector similar in technology to 
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Intime cosmics/yr Depth (mwe) 



5x 10^ 





4230 


1050 


462 


2000 


77 


3000 


15 


4400 



TABLE IX: Number of cosmic ray muons in a 50 m height/diameter detector in a lO^s pulse for 10^ pulses, 
corresponding to approximately 1 year of running, versus depth in meters water equivalent. 



Super- Kamiokande, the dead-time from the above event rates will exceed 50% [|35l . To reduce this 
dead-time using fast pulse digitizers is costly, and requires significant software and hardware R&D 
to resolve overlapping pulses to reconstruct events with contained vertices. The consequences on 
background rejection and resolution are at present unknown. The depth required to reduce the 



number of in-time cosmics to various levels is given in Table IX A depth of at least ~ 1000 meters 
water equivalent is needed to reduce the muon rate to a level comparable to the rate of events from 
the neutrino beam so that minimal dependence on pattern recognition (and a modest active veto 
capability) is needed to separate beam related events. 

A 50 kT liquid argon TPC can be contained in a cylindrical tank of size 35.5 m height/diameter; 
such a detector will have a cosmic ray muon rate of 125 kHz from the top and 125 kHz from the 
sides. An examination of cosmic rays [iTOl in a liquid argon TPC has considered their effects on 
data acquisition and event reconstruction, and as a source of background. The rate of cosmic rays 
was shown to be tolerable with the proposed drift-time and data acquisition system for cycles up 
to 5 Hz. In this scheme the detector takes data in a short time interval (currently proposed to be 
3 drift times) near the beam time. This is sufficient to cover most possible accelerator cycle times 
discussed above. The high granularity of the detector should allow removal of cosmic muons from 
the data introducing a small (< 0.1%) inefficiency to the active detector volume, so that most of 
the accelerator induced events are unobscured. If a cosmic ray muon (photon) event mimics a 
contained in-time neutrino event it must be rejected based on pattern recognition. The rejection 
required is ~ 10^ for muon cosmics and ~ 10^ — 10^ for photon cosmics; given the fine grained 
nature of the detector this rejection is considered achievable, but still needs to be demonstrated by 
detailed simulations. 
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13. ANALYSIS OF SENSITIVITY TO OSCILLATION PARAMETERS 



In this section we will combine the information from the previous sections on the intensity of 
the accelerator beam and detector performance to calculate the sensitivity to oscillation parameters. 
The main features of the accelerator and detector performance can be summarized as follows: 

For the sensitivity calculations we will assume that we can obtain a total of 60 x 10^^ protons 
at 120 GeV. This total is to be divided between neutrino and anti-neutrino running. To convert 
this luminosity to running time we will assume that the accelerator can produce proton intensity 
according to Figure |7] in the accumulator upgrade scenario. In the accumulator upgrade scenario 
a power level of 1.2 MW is expected at 120 GeV. This running scenario will be used for both 
the off-axis and the DUSEL based options for the Phase II program. We will make comments on 
running at lower proton energies as well as more exposure. The impact on the running time will 
be according to the power curve in Figure |7J The raw event rates can be obtained from Tables IV 
tolVIIl 

For the first DUSEL based calculation we have assumed a water Cherenkov detector with a 



total fiducial mass of 300 kTon with the performance described in Section 10.1 The calculation 
was performed with the GLoBES package [44] with the beam spectra and detector performance 
specified according to the work in this report. For the DUSEL baseline we have performed cal- 
culations ranging from 500 km to 2500 km with various beam configurations. We cannot display 
aU calculations in this report due to length considerations, but they can be obtained from the study 
website[|15l. Differences in parameter sensitivity due to baseline will be discussed. Most of the 
calculations shown here will be for the 1300 km distance. 

For the second DUSEL based calculation we have assumed a 100 kTon liquid argon time pro- 



jection chamber with the performance indicated in section 10.2 Briefly, we assume 80% efficiency 
for electron neutrino events with very little background from other sources. 

For the off-axis calculations several different combinations were calculated. First, for com- 
parison purposes the calculation is performed for NOvA with the detector performance obtained 
from the NOvA collaboration. Second, a 100 kTon total mass for a liquid argon detector TPC was 
assumed for phase II. The performance was evaluated for setting the entire detector mass at the 
same location as NOvA and also for setting 50 kT at the NOvA site (12 km off-axis) and 50 kT at 
the site (40 km off-axis) where the second oscillation maximum can be observed. Doubling of the 
total mass at the two sites was also examined. 

Lastly, we note that unless otherwise noted the oscillation parameters used for the calculations 
are as follows: 

Am^i = 8.6 X IQ-^eV^ 
sin^ 2012 = 0.86 
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Am^2 = 2.7 X lO^^eV^ 
sin2 2023 = 1 

The parameters, 613, dcp and the mass hierarchy (normal or reversed) are left free in the calcu- 
lation. 

Before describing the sensitivity to Vg appearance, we first make a few comments on the 
disappearance measurement. Either of the two experimental concepts (a new beam to a DUSEL 
location or new off axis detectors on the surface in the NuMI beam-line) have sufficient statistical 
reach to make a very precise measurement of the atmospheric oscillation parameters (Amgj and 
sin 2023- For a 100 kTon detector, ~ 10000 CC events per year are expected in either sce- 
nario with approximately 1/2 disappearing due to oscillations (see Section [8]for exact numbers for 
specific beam configurations). The statistical precision after several years of running, therefore, 
will be < 1% for both Am32 and sin 2023- A discussion of this measurement for a DUSEL based 
detector can be seen in [|20ll . A similar discussion for the off-axis scenario is in 11721 . To obtain the 
best measurement, both the neutrino event energy resolution (including nuclear target effects due 
to Fermi motion and re- scattering) and the absolute energy scale need to be well modeled. With 
current knowledge of these limitations the measurement will most likely be systematically limited 
to about 1% for both Am32 and sin^2023- 

There is an important difference between the off-axis measurement and the broad band mea- 
surement. The oscillation shape including a nodal pattern, if the baseline distance is sufficient, can 
be measured with a DUSEL based detector. Such a measurement will exhibit less correlation be- 
tween the two parameters Am32, which determines the position of the node in energy, and sin 2023, 
which determines the depth (or amplitude) of the node. A precise measurement of the shape could 
also limit non-standard physics models of decay, decoherence, extra-dimensions, etc. 

13.1. Sensitivity of a FNAL to DUSEL based program 

For the calculations reported here we have used the 120 GeV beam with 380 meter decay tunnel 
with a 0.5 deg off-axis angle. As explained above, the spectrum from such a configuration is well 
matched to the physics at this current time. The energy of the proton beam and the horn optics 
need to be optimized further. 
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13. 1. 1. Water Cherenkov Detector 



The reconstructed electron neutrino spectrum with 300 kT of fiducial mass and a total exposure 
of 60 X 10^*^ protons (divided equally between neutrinos and antineutrino running) is shown in 
Figure |9} This spectrum includes effects of nuclear motion, detector resolution, detector signal 



efficiency, and background rejection using the performance as described in Section 10.1 The plot 
is made for sin 20i3 = 0.04 and a baseline of 1300 km. The plots for 1480 km can be obtained 
from [[151. The left plots are for the normal mass hierarchy (nii < m2 < m^). The right hand plots 
are for the reversed hierarchy (itit, <m\ < ml). The top plots are for neutrino running and bottom 
plots are for anti-neutrino running. 

By fitting the spectra in Figure|9]we can extract the parameters 613, 5cp, and the mass hierarchy. 
We calculate a function and extract the confidence levels for a simultaneous fit to these three 
parameters. For the input values of the other oscillation parameters we assume 1 sigma errors as 
follows: 

0i2 = O.59±lO%,Am^i = (0.86 ±10%) x 10"^ 

023 = ;r/4±5%,Am^i = (2.7 ±5%) x 10"^ 

We also include 5% error on the matter density. The calculation includes correlations between 
all parameters and accounts for possible degeneracies. The spectra were fit with statistical errors 
and with 10% systematic error on the background and 1% systematic error on the normalization 
with no correlations between neutrino and anti-neutrino channels. Details of the analysis method 
are in [[731 where the same analysis was performed with a different spectrum and detector perfor- 
mance. 



In Figure 10 we show the confidence level contours for measuring the pair of parameters (613 
and 5cp)- This calculation was performed for normal mass hierarchy, a baseline of 1300 km, and 
a total exposure of 60 x 10^^ protons equally divided between neutrino and anti-neutrino running. 
The result for 1480 km is approximately the same. In the case of normal hierarchy, the neutrino 
data alone can be used to measure the parameters over a large range of parameter space. But if 
the mass hierarchy is reversed, anti-neutrino data has to be used. The resolution obtained after 
combining both neutrino and anti-neutrino data is approximately independent of mass hierarchy. 
It is clear that the parameter measurement will suffer from background below sin 20i3 = 0.01, 
but above this value the resolution on the CP phase of about ±20"^ (1 sigma) is approximately 
independent of ^13. 

'-) 

If there is no excess of electron events observed then we can set a limit on the value of sin 2 613 



as a function of 5cp- Such sensitivity limits are shown in Figure 1 1 The range of parameters over 
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which the mass hierarchy can be resolved is shown in Figure 12 We have chosen to display the 
limits separately for the two mass hierarchies. Some of the structure in the 3 sigma lines is due 
to the hmited number of bins used in the calculation. The region to the right hand side of each 
curve excludes the opposite mass hierarchy at the respective confidence level. Similarly the range 
of parameters over which CP violation can be established (i.e. determine that 5cp is not or n) is 



displayed in Figure 13 



53 



> : V, 30 1(f°PoT, 1300km 
^ normal hierarchy 
Ln - sin^2e,, = 0.04 



signal + bkg: 



5(,p=+45° (702) 
8cp=0° (807) 
45° (934) 
background: 

all (415) 
beamVg (196) 




neutrino energy [GeV] 



> : V, 30 ltf°PoT, 1300km 

45- normal hierarchy 
Ln : sin^2e., = 0.04 



signal + bkg: 
— 5^p=+45° (366) 
♦ 8,p=tf> (342) 
45° (311) 
background: 

all (201) 
beamVg (121) 




> 
a> 
(D 
m 

(M 

o 

> 
LU 



90 
80 
70 
60 
50 
40 
30 
20 
10 


50 



V, 30 ^&°Po^, 1300km 
_ inverted hierarchy 
sin^2e, =0.04 



signal + bkg: 

— 5^p=+45° (538) 
t S,p=0° (607) 

- S^p=-45° (688) 
background: 

all (419) 
beamVg (199) 




neutrino energy [GeV] 



O 45 



10 

(M 



40 



h V, 30 itf^PoT, 1300km 

inverted hierarchy 
t sin^ 26^3 = 0.04 



t 35 



LU 



30 
25 
20 
15 
10 



signal + bkg: 
- S,p=+45°(493) 

8cP=0" 

45° (395) 

background: 

all (200) 
beamVg (119) 




neutrino energy [GeV] 



neutrino energy [GeV] 



FIG. 9: Simulation of detected electron neutrino (top plots) and anti-neutrino (bottom plots) spectrum (left 
for normal hierarchy, right for reversed hierarchy) for 3 values of the CP parameter 5cp, —45°, 0", and 
—45", including background contamination. This simulation is for 300 kT of water Cherenkov detector with 



the performance described in Section\10.1\ This is for an exposure of 30 x 10^" POT for each neutrino and 



anti-neutrino running. The hatched histogram shows the total background. The beam background is also 
shown. The other parameters and running conditions are shown in the figure. 
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FIG. 10: 90% and 95% confidence level error contours in sin^ 20i3 versus dcpfor statistical and systematic 
errors (left hand plot) for 15 test points. This is for a 300 kT water Cherenkov detector with a total exposure 
of 60 X lO^'' POT. The right hand side is for statistical errors alone. This plot was made for normal mass 
hierarchy. We assume 10% systematic errors on the background for this plot. 
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FIG. 11: 3 sigma and 5 sigma confidence level exclusion limits for determining a non-zero value for 613 
in sin^20i3 versus dcp for statistical and systematic errors (left hand plot). This is for a 300 kT water 
Cherenkov detector with a total exposure of 60 x 10^" POT. The right hand side is for statistical errors 
alone. The solid (dashed) lines are for normal (reversed) mass ordering. We assume 10% systematic errors 
on the background for this plot. 
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FIG. 12: 3 sigma and 5 sigma confidence level exclusion limits for determining the mass hierarchy in 
sin^26i3 versus 5cp for statistical and systematic errors (left hand plot). This is for a 300 kT water 
Cherenkov detector with a total exposure of 60 x 10^" POT. The right hand side is for statistical errors 
alone. The solid (dashed) lines are for normal (reversed) mass ordering. We assume 10% systematic errors 
on the background for this plot. 




sin^2e,g sin^2e^3 

FIG. 13: 3 sigma and 5 sigma confidence level exclusion limits for determining CP violation in sin^20i3 
versus dcpfor statistical and systematic errors (left hand plot). This is for a 300 kT water Cherenkov detector 
with a total exposure of 60 x lO^'' POT. The right hand side is for statistical errors alone. The solid ( dashed) 
lines are for normal ( reversed) mass ordering. We assume 10% systematic errors on the background for this 
plot. 
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Sensitivity variation with exposure: The exposure assumed in the above plots was 30 x 10^*^ 
protons on target for each neutrino and antineutrino running. This corresponds to 3 years of running 
for each polarity for 1.2 MW of beam power and 1.7 x 10^ sec per year of running at 120 GeV. 
If we were to run the antineutrino beam for twice the exposure of neutrino, then the the number 
of events is approximately balanced between v and v. Such unequal running is advantageous in 
the case of the reversed hierarchy. An analysis of the total exposure was performed in [|73ll . It was 
found that longer exposures will have relatively modest effect on the sensitivity to sin 20i3 and the 
mass hierarchy resolution, but could be important for improving the precision on the CP violation 
measurement. Exclusion contours for twice the exposure (total exposure of 120 x 10^*^ protons) 
are shown in Figures [14] (for determining non-zero 613), 15 (for determining mass hierarchy), and 
T6](for determining CP violation). 
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FIG. 14: i sigma and 5 sigma confidence level exclusion limits for determining a non-zero value for 613 
in sin^20i3 versus dcp for statistical and systematic errors (left hand plot). This is for a 300 kT water 
Cherenkov detector with a total exposure of 120 x lO^'' POT. The right hand side is for statistical errors 
alone. The solid (dashed) lines are for normal (reversed) mass ordering. We assume 10% systematic errors 
on the background for this plot. 

Sensitivity variation with distance: Analysis in fTSl showed that there is significant variation 
in sensitivity up to 1500 km for determination of the mass hierarchy. This is reproduced in Figure 
17 The variation in sensitivity to 613 was found to be mild partly because of the larger matter 
enhancements (in neutrino (antineutrino) mode for normal (reversed) mass hierarchy) at longer 
distances. There is a slow decrease in the sensitivity to CP violation at longer distances, but this is 
attributed to the shape of the spectrum used for this calculation. At distances below 1000 km, there 
is a degradation in the CP sensitivity because of the need to resolve the ambiguity due to the mass 
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FIG. 15: 3 sigma and 5 sigma confidence level exclusion limits for determining the mass hierarchy in 
sin^26i3 versus 5cp for statistical and systematic errors (left hand plot). This is for a 300 kT water 
Cherenkov detector with a total exposure of 120 x lO^'' POT. The right hand side is for statistical errors 
alone. The solid (dashed) lines are for normal (reversed) mass ordering. We assume 10% systematic errors 
on the background for this plot. 



hierarchy with the same data. Complete calculations with same assumptions for detector size and 
performance and spectra can be obtained from the [fTSl for distances up to 2600 km; this covers 
the various options for the DUSEL locations. 

Sensitivity variation due to systematics and parameter variation: 

The sensitivity calculation reported in this section follows the prescription from [14411731 . They 



include the parameter variation as described above (Section 13.1.1 1. In addition, we assume a 10% 
systematic error on the total background. Considering recent and past experience with background 
determination in long baseline experiments, the 10% systematic error is very likely a pessimistic 
assumption [|74l 1751 . especially with a planned near detector. Sensitivity estimates with other as- 
sumptions for the systematic error can be obtained from the website [fTSl . The conclusion from 
these studies is that the background systematic error will most likely dominate over the parameter 
variation. Therefore, in Figures [10] to [16} we have chosen to show the sensitivity with and without 
systematic errors. Our conclusion is that the wide band technique which leads to a spectrum mea- 
surement is robust against parameter changes and background systematic errors over a reasonable 
range fTSl . 
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FIG. 16: 3 sigma and 5 sigma confidence level exclusion limits for determining CP violation in sin^20i3 
versus 5cp for statistical and systematic errors (left hand plot). This is for a 300 kT water Cherenkov 
detector with a total exposure of 120 x 10^" POT. The right hand side is for statistical errors alone. The 
solid (dashed) lines are for normal (reversed) mass ordering. We assume 10% systematic errors on the 
background for this plot. 
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FIG. 17: Discovery reach for a normal mass hierarchy at 3a for CP fractions (lower-most line, best case), 
0.5 (middle line) and 1 (uppermost line, worst case) as a function of the baseline. The detector mass, beam 
power and exposure are kept the same for all baseUnes. For further explanation of the plot please see fTBll . 



59 



13. 1.2. Liquid Argon Detector 



If a 100 kTon fiducial mass liquid argon time projection chamber can be built, then it can 
be placed at one of the DUSEL sites and used as a long baseline neutrino oscillation detector. 
We have assumed that such a detector can have 80% efficiency for all charged current electron 



neutrino events and has background rejection capability (Section 10.2 1 that virtually rejects all NC 
and CC backgrounds. We further assume that the detector will have resolution characterized by 



20% / a/ GeV for non-quasielastic events and 5%/ a/^/ GeV for quasielastic events. The spectra 
that result from these assumptions are displayed in Figure[T8]for the same parameters and exposure 
as Figure |9j 

The parameter resolutions and sensitivity limits for the 100 kT liquid Argon TPC at DUSEL 



are shown in Figures 19 to 22 If there is no excess of electron events observed then we can set 
a limit on the value of sin^20i3 as a function of 5cp. Such sensitivity limits are shown in Figure 



20 The range of parameters over which the mass hierarchy can be resolved is shown in Figure 



21 We have chosen to display the limits separately for the two mass hierarchies. The region to 



the right hand side of each curve excludes the opposite mass hierarchy at the respective confidence 
level. Similarly the range of parameters over which CP violation can be established (i.e. determine 



that 5cp is not or it) is displayed in Figure 22 The comments we made regarding dependence 
on exposure, baseline, and oscillation parameters for the water Cherenkov detector are equally 
applicable to the sensitivities one would obtain with the liquid argon detector placed at DUSEL. 
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FIG. 18: Simulation of detected electron neutrino (top plots) and anti-neutrino (bottom plots) spectrum 
(left for normal hierarchy, right for reversed hierarchy) for 3 values of the CP parameter 5cp, —45°, 0", 
and —45", including background contamination. This simulation is for 100 IcT of LAr detector (with the 
performance described in the text) placed at DUSEL, 1300 km away from FNAL. The hatched histogram 
shows the total background, which is dominated by the Vg beam background. The other parameters and 
running conditions are shown in the figure. 
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FIG. 19: 90% and 95% confidence level error contours in sin^ 20i3 versus dcpfor statistical and systematic 
errors (left hand plot) for 15 test points. This plot is for a 100 kTon liquid Argon TPC placed at DUSEL 
1300 km away from FNAL. This is for combining both neutrino and anti-neutrino data. The right hand side 
is for statistical errors alone. This plot was made for normal mass hierarchy. We assume 10% systematic 
errors on the background for this plot. 




sin^2e^3 sin^2e^3 



FIG. 20: 3 sigma and 5 sigma confidence level exclusion limits for determining a non-zero value for 0i3 in 
sin^20i3 versus 8cp for statistical and systematic errors (left hand plot). This plot is for a 100 kTon liquid 
Argon TPC placed at DUSEL 1300 km away from FNAL. This is for combining both neutrino and anti- 
neutrino data. The right hand side is for statistical errors alone. The solid (dashed) lines are for normal 
( reversed) mass ordering. We assume 10% systematic errors on the background for this plot. 
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FIG. 21: 3 sigma and 5 sigma confidence level exclusion limits for determining the mass hierarchy in 
sin^26i3 versus 5cp for statistical and systematic errors (left hand plot). This plot is for a 100 kTon liquid 
Argon TPC placed at DUSEL 1300 km away from FNAL. This is for combining both neutrino and anti- 
neutrino data. The right hand side is for statistical errors alone. The solid (dashed) lines are for normal 
(reversed) mass ordering. We assume 10% systematic errors on the background for this plot. 
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FIG. 22: 3 sigma and 5 sigma confidence level exclusion limits for determining CP violation in sin^20i3 
versus dcpfor statistical and systematic errors (left hand plot). This plot is for a 100 kTon liquid Argon TPC 
placed at DUSEL 1300 km away from FNAL. This is for combining both neutrino and anti-neutrino data. 
The right hand side is for statistical errors alone. The solid (dashed) lines are for normal (reversed) mass 
ordering. We assume 10% systematic errors on the background for this plot. 
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13.2. Sensitivity of a NuMI based off axis program 



In this section we present the sensitivity results for several of the scenarios using NuMI off-axis 
neutrinos which we studied. The following assumptions were made for this part of the study. We 
have attempted to make the calculations as directly comparable to the previous section as possible. 
But there are small differences. These are also listed below. 

• The source of the proton beam used to create neutrinos will be the Fermilab Main Injector. 

• Current planning at Fermilab to maximize the proton intensity from the Main Injector via a 
series of staged upgrades, will in fact occur over the next decade. 

• The possibility of a new source, such as the HINS (High Intensity Neutrino Source [,25 J ), or 
alternative ideas to improve the Main Injector, is considered as the final stage of the upgrade 
path, ultimately providing an annual proton intensity of 2 x 10^^ protons per year. 

• Though the Fermilab Main Injector can be operated with extracted proton energies of less 
than 120 GeV, we have assumed that the optimum operation is at 120 GeV, based on total 
delivered beam power. 

• The neutrino beam which we are considering here is the existing NuMI beam, which is 
a conventional horn focusing beam, capable of producing both neutrino and anti-neutrino 
beams (by reversing the current in the horns). We do not consider any reconfiguration or 
modification to the existing 2-meter diameter, 675-meter long decay pipe. 

• We assume that upgrades to targets, horns, shielding and cooling systems will be required 
to accommodate proton intensities significantly higher than that for which the facility was 
designed. 

• We do assume that the target and horn configuration can be adjusted to optimize neutrino 
rates. 

• We assume that detectors situated in NuMI off-axis locations will most likely be sited on or 
near the surface. We do not discuss the detector designs required to reject the backgrounds 
from cosmic ray interactions associated with a surface location. 

• We have used an efficiency of 80% for all charged current electron neutrino events and a 
neutral current rejection factor of 0.001, consistent with the parameters of a Liquid Argon 
detector. 

• We assume that the background is known with a systematic error of 5%. 
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• The assumption on oscillation parameters is stated at the beginning of this section. For the 
calculation here these parameters as well as the matter density are assumed to be fixed. 

• In our first pass analysis our sensitivities have been generated assuming no prior knowledge 
of either sin^20i3 , the mass hierarchy or dcp ■ Iteration of the sensitivity calculation for 
the mass hierarchy and dcp have also been done such that for values of sin^20i3 > 0.02, the 
angle is known (as will be the case from Phase I experiments). For values < 0.02 we assume 
the angle is unknown. 

As a cross-check and starting point for our study, we have calculated the sensitivities for a 20 



kton NOvA detector (see Figure 23 ). Our results are consistent with those produced for the NOvA 
project Technical Design Report. Note, for this and all subsequent figures, the dashed line is placed 
at the current Chooz limit. 

For each scenario we assumed that the Phase II program consisted of running for an equal 
time in neutrino and anti-neutrino mode. We show the plots for an integrated proton intensity of 
30 X 10^*^ in each mode. We also assume that the NOvA detector continues to take data during 
Phase II. The new detectors are all Liquid Argon technology. 

The, first scenario considered was placing a 100 kton detector at the 1st maximum, i.e. simply 
increasing the mass and efficiency of the NOvA configuration. The sensitivities are shown in 



Figure 24 



The second scenario we studied, was to place a 100 kton detector at a baseline of 700 km and an 
off-axis angle of 57 mrad (40 km). This location corresponds to the second oscillation maximum, 
where the matter effects are small (due to the lower energy of the neutrinos), but the CP effects 
are large. The NOvA detector is the only detector at the first maximum site (L = 810, 14 mrad off 



axis). These results are shown in Figure 25 We find that with this configuration, running neutrinos 
and anti-neutrinos, the sensitivity to the mass hierarchy flattens over the range of possible dcp 
values, but the discovery potential is limited to values of 613 relatively close to the current limit. 

A third scenario was to split the mass between the two locations, 50 ktons at each the first and 
second maximum. These results are in Figure [26| 

Finally, di fourth scenario was the same as the third except that the detectors at each site were 
100 ktons. 

We summarize our studies in Table 1X1 Because some of the scenarios studied have the benefit 
of "flattening" the sensitivity over 5cp , we have included the sensitivity limits for both 50% and 
100% coverage of the 5cp space. We have also included in this table the limits which can be 
reached as the Phase I program evolves, and as the Liquid Argon technology also evolves. 

From these studies we conclude the following: 
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• sin^20i3 down to 0.02 can be measured by the Phase I (NOvA) experiment. Phase I ex- 
periments however, have limited or no sensitivity to determining the mass hierarchy, and 
essentially no sensitivity to 5cp ■ 

• If sin^20i3 is large, i.e. >0.04 a Phase II experiment using the NuMI beam can be designed 
specifically to determine the mass hierarchy. Such an experiment is like our 4th scenario 
(see Table [X]), with two massive detectors placed at the first and second maximum sites. In 
this experiment the mass hierarchy can be resolved for all values of 5cp ■ The experiment 
also has sensitivity to 5cp . 

• The most interesting and complex situation to plan for, is if Phase I experiments indicate that 
0.02 < sin^20i3 < 0.04. In this case we find that the configurations studied for the NuMI 
Off- Axis option can have relatively good sensitivity to determine the mass hierarchy, as well 
as some sensitivity to CP. 

• If Phase I experiments conclude that sin^20i3 < 0.02 the Phase II program can continue the 
search. Continued running, more protons and larger more efficient detectors, placed at the 
1st maximum (the NOvA site), allows one to reach sensitivities to well below 10^^, (of the 
order ~0.003) as can be seen in our first scenario. 



13.3. Comparison of sensitivity estimates 

A summary of the sensitivity reach for non-zero 613, CP violation and the sign of Amgj for 
6 different combinations of beams, baselines, detector technologies, and exposure is presented in 



Table |XIj Several more configurations for the off-axis scenario are presented in Table |Xj The 
sensitivity reach is given as the lowest sin^ 20i3 value at which at least 50% of 5cp values will have 



> 3a reach. For this table we use the mass hierarchy with the worst sensitivity to determine the 
minimal value of sin^ 2 613 for which > 50% of 5cp values will have > 3c7 sensitivity to a particular 

'-) 

measurement. We estimated these values of sin 2 613 from the studies and plots discussed in 



Sections 13.1 and 13.2 We note that different options are sensitive to different values of 8cp, such 
that being sensitive to 50% 5cp values does not necessarily imply that a given experimental option 
is sensitive to the same region of oscillation parameter phase space as another. 

We compare the wide-band FNAL to DUSEL program, option (4), with the narrow-band off- 
axis NuMI-based program, option (2), for the same exposure of 6.8 MW.yr (1 experimental year is 
defined as 1.7 x 10^ seconds). This is equivalent to an integrated exposure of 60 x 10^*^ protons-on- 
target for proton beam energies of 120 GeV. We assume equal amounts of exposure for neutrinos 
and anti-neutrino (reverse horn current) running. A liquid Argon TPC with a total mass of 100 
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kT is assumed as the detector technology of choice for the purpose of the comparison. We note 
that slightly different assumptions on the systematic uncertainties on the oscillation parameters and 
backgrounds went into the sensitivity estimates for NuMI off-axis (5% uncertainty on the back- 
ground) and the wide-band FNAL to DUSEL options (10% uncertainty on the background). The 
effect of the different assumptions is < 15% variation on the value of sin^ 20i3 at which the sensi- 
tivity reaches 50% of dcp- We find that for the same exposure of 6.8 MW.yr, and the same liquid 
Argon TPC detector technology (size and same performance), the wide-band FNAL to DUSEL 
approach has significantly better sensitivity to CP violation, the sign of Amg^ and comparable 
sensitivity to non-zero values of 613. To illustrate the improvement in sensitivity over the existing 



program, the sensitivities of the current NOvA experiment (as shown in Figure 23 1 at the same 



exposure, are summarized as option (1) in Table XI These are the sensitivity limits expected from 
the NOvA experiment only, before combination with T2K ir76l . Analysis of combining with T2K 
has been performed elsewhere iTTTll . 

The value of sin^ 2 613 at which at least 50% of 5cp values will have > reach as a function of 
exposure for the NuMI ME beam at 810km (labeled NOvA*), and the wide-band 120 GeV beam 
at 1300km (labeled WBB-120s) is summarized in Figure |28]from reference [|78l . A LAr TPC 
is the detector technology assumed for NuMI off axis (labeled NOvA*) and WBB-120s. We find 
that after reaching an exposure of 2 MT.MW.IO^ seconds (for lOOkT LAr and a 120 GeV beam, 
this is an exposure of 10^^ protons -on-target), the mass hierarchy-5cp degeneracy is sufficiently 
resolved for the NOvA* approach (option (2)) - and the sensitivity to CP violation approaches that 
of the wide-band beam at the 1300km baseline. For the mass hierarchy, the wide -band FNAL to 
DUSEL approach always has significantly better sensitivity independent of the exposure. Option 



(3) in Table XI is a NuMI-based program with a 50 kT detector at the 1st oscillation maximum 
running concurrently with another 50kT module placed at the 2nd oscillation maxima. We find 
that option (3) has worse sensitivity to non-zero values of 613 when compared to option (2) and 
slightly better sensitivity to the sign of Awgj. 

Option (5) summarizes the FNAL to DUSEL sensitivity when the 100 kT LAr TPC of option (4) 
is replaced by a 300 kT water Cherenkov detector. We find that the sensitivity worsens due to the 
lower signal statistics and higher NC backgrounds in a water Cherenkov detector. We can recover 
some of the lost sensitivity by doubling the exposure of the water Cherenkov detector as shown in 
option (6). For the same exposure, the FNAL to DUSEL program with a 300 kT water Cherenkov 
detector, option (5), has the same sensitivity to CP violation as the NuMI based program with a 100 
kT LAr TPC in options (2) and (3) and significantly better sensitivity to the sign of Amjj. We find 
the FNAL to DUSEL program with a 300 kT water Cherenkov detector has similar sensitivity to 
non-zero 613 as the NuMI based program with two 50 kT LAr TPC's at the 1st and 2nd oscillation 
maxima, option (3). 
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We summarize the comparison studies as follows: 

• sin^20i3 down to 0.02 can be measured by the Phase I (NOvA) experiments. Phase I 
(NOvA) experiment, however, have limited or no sensitivity to determining the mass hi- 
erarchy, and essentially no sensitivity to 5cp . 

• All Phase II experimental options will improve the sensitivity to CP violation by at least an 
order of magnitude over the existing Phase I program. 

• Given the same exposure and detector technology (LAr TPC), the FNAL to DUSEL program 
with a wide band beam has significantly better overall sensitivity to neutrino oscillations 
when compared to a shorter baseline NuMI based program with an off-axis beam. 

• The FNAL to DUSEL program with a 300 kT water Cherenkov detector has similar sensi- 
tivity to CP violation when compared to a NuMI off-axis program with a 100 kT LAr TPC, 
and significantly better sensitivity to the sign of Amjj. 

• A NuMI off-axis program with two 50 kT LAr TPCs at the 1st and 2nd oscillation maxima 
at baselines of 810 and 700 km respectively has marginally better sensitivity to the sign of 
Arn^i but significantly worse sensitivity to non-zero O12 when compared with putting the full 
100 kT mass at the 1st oscillation maxima. 



68 




FIG. 23: 90%, 3 sigma confidence level exclusion limits for determining a non-zero value for (top), for 
excluding the opposite mass hierarchy (middle), and for excluding CP violation (bottom) in sin^20i3 versus 
dcp- These plots (blue for normal and red for reversed hierarchy) are for a 20 IcTon NOvA detector placed 
at the off-axis location on the NuMI beam-line with a total exposure of 60 x lO^'' protons, and for combining 
both neutrino and anti-neutrino data. 5% background systematic errors are assumed. 
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FIG. 24: Scenario 1 : 90%, 3a, and 5a confidence level exclusion limits for determining a non-zero value 
of 013 (left), for excluding CP violation (center), and for excluding the opposite mass hierarchy (right), in 
sin^20i3 versus dcp ■ 




FIG. 25: Scenario 2 : 90%, 3a, and 5a confidence level exclusion limits for determining a non-zero value 
of 013 (left), for excluding CP violation (center), and for excluding the opposite mass hierarchy (right), in 
sin^20i3 versus 5cp ■ 
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FIG. 26: Scenario 3 : 90%, 3a, and 5a confidence level exclusion limits for determining a non-zero value 
of 013 (left), for excluding CP violation (center), and for excluding the opposite mass hierarchy (right), in 
sin^20i3 versus dcp ■ 




FIG. 27: Scenario 4 : 90%, 3a, and 5a confidence level exclusion limits for determining a non-zero value 
of 013 (left), for excluding CP violation (center), and for excluding the opposite mass hierarchy (right), in 
sin^20i3 versus 5cp ■ 
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Detector Location 
(L(km), 0(mr)) 
(beam tune) 


Detector 
(technology/mass) 


Exposure 
POT(V) 
/POT(V) 


sin220i3 


sgn(Amf 1 ) 


CPV 


810, 14(ME) 


NOvA/20kt 


15/- 


0.018/0.030 


0.17/NA 


NA 


810, 14 (ME) 


NOvA /20kt 


15/15 


0.018/0.024 


0.16/NA 


NA 


810, 14(ME) 


NOvA /20kt 


30/30 


0.012/0.020 


0.10/NA 


NA 


810, 14(ME) 
810,14(ME) 


NOvA /20kt 
+ LAr/5kt 


30/30 
30/30 


0.007/0.013 


0.08/0.20 


NA 


810, 14(ME) 
810,14(ME) 


NOvA /20kt 
+ LAr/20kt 


30/30 
30/30 


0.004/0.009 


0.05/0.15 


0.07/NA 


Scenario 1 
810, 14(ME) 
810, 14(ME) 


NOvA/20kt 
+ LAr/lOOkt 


30/30 
30/30 


0.0018/0.005 


0.03/0.12 


0.03/NA 


Scenario 2 
810, 14(ME/LE) 
700, 57 (LE) 


NOvA/20kt 
+ LAr/lOOkt 


30/30 
30/30 


0.011/0.018 


0.05/0.07 


0.07/NA 


Scenario 3 
810, 14(ME/LE) 
810, 14 (LE) 
700, 57 (LE) 


NOvA/20kt 
+ LAr/ 50kt 
+ LAr/ 50kt 


30/30 
30/30 
30/30 


0.0035/0.006 


0.033/0.06 


0.035/NA 


Scenario 4 
810, 14(ME/LE) 
810, 14(LE) 
700, 57 (LE) 


NOvA /20kt 
+ LAr/lOOkt 
+ LAr/lOOkt 


30/30 
30/30 
30/30 


0.0027/0.0046 


0.030/0.042 


0.022/NA 



TABLE X: Sensitivity comparisons for all NuMI Off-axis scenarios that were evaluated. These numbers 
were calculated with the normal hierarchy assumption. The first three cases represent three possible stages 
of the Phase I (NOvA ) program. The values given represent the value of sin^20i3 where a 3 C7 determination 
of the parameter can be made for 50%(/ 100%) of the possible values of 5cp . Note that for determining the 
sensitivity to mass hierarchy and 5cp , for values of sin^20i3 >0.02 we assume that the angle is known (i.e. 
from Phase I experiments.) 
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Option 


Beam 


Baseline 


Detector 


Exposure (MW.yr*) 


013/0 


CPV 


/A 9 \ 

sgn{Am^y) 


(1) 


NuMI ME, 0.9° 


810 km 


NOvA 20 kT 


6.8 


0.015 


>0.2 


0.15 


(2) 


NuMI ME, 0.9° 


810 km 


LAr 100 kT 


6.8 


0.002 


0.03 


0.05 


(3) 


NuMI LE, 0.9°, 3.3°, 


810,700 km 


LAr 2 X 50 kT 


6.8 


0.005 


0.04 


0.04 


(4) 


WBLE 120GeV, 0.5° 


1300km 


LAr 100 kT 


6.8 


0.0025 


0.005 


0.006 


(5) 


WBLE 120GeV, 0.5° 


1300km 


WCe 300 kT 


6.8 


0.006 


0.03 


0.011 


(6) 


WBLE 120GeV, 0.5° 


1300km 


WCe 300 kT 


13.6 


0.004 


0.012 


0.008 



TABLE XL Comparison of the sensitivity reach of different long baseline experiments. The sensitivity is 
given as the value of sin^20i3 at which 50% of dcp values will have > 3a reach for the choice of mass 
hierarchy with worst sensitivity. We assume equal amounts of v and v running in the total exposure. The 
assumption on running time is 1.7 x 10^ seconds of running per year. Also see Table [x| 
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sin^2ei3 



NOvA*- 
WBB-120s' 
T2KK- 




FIG. 28: The sin^20i3 reach at 3a for the discovery of nonzero sin^20i3 (top), CP violation(middle), and 
the normal hierarchy (bottom) as a function of exposure. The curves are for a fraction of 5cp of 0.5, which 
means that the performance will be better for 50% of all values of 5cp, and worse for the other 50%. The 
light curves in the CPV panel are made under the assumption that the mass hierarchy is known to be normal. 
The shaded regions result by varying the systematic uncertainties from 2% (lower edge) to 10% (upper 
edge). This figure is reproduced from fTSll . 
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14. SENSITIVITY TO NON-ACCELERATOR PHYSICS 



A well instrumented very large detector, in addition to its accelerator based neutrino program, 
could be sensitive to nucleon decay which is one of the top priorities in fundamental science. All 
of the detector technologies we consider will lead to enhanced detection and study of neutrinos 
from natural sources such as the Sun, Earth's atmosphere and lithosphere, and past and current 
supernova explosions. To achieve these goals, the key issues are cosmogenic backgrounds and 
low energy thresholds (~ 5 MeV); the first primarily depends on the depth of the detector and 
the second depends on the depth, the radioactivity from the materials used in the detector and in 
the surrounding rock, and the detector noise (photosensor noise in the case of a water Cherenkov 
detector, and electronic noise in the case of a liquid argon TPC). 

In this section we briefly summarize the potential of a large detector for nucleon decay and 
astrophysical sources of neutrinos. We also comment on the technical requirements on the detector. 
For each topic we attempt to identify where the requirements of the accelerator program match and 
where they diverge. 

14.1. Improved Search for Nucleon Decay 

Theoretical Motivation: While current experiments show that the proton lifetime exceeds about 
j^q33 ygars^ its ultimate stability has been questioned since the early 1970's in the context of theo- 
retical attempts to arrive at a unified picture of the fundamental particles - the quarks and leptons 
- and of their three forces: the strong, electromagnetic and weak. These attempts of unification, 
commonly referred to as "Grand Unification", have turned out to be supported empirically by the 
dramatic meeting of the strengths of the three forces that is found to occur at high energies in 
the context of so-called "Supersymmetry", as well as by the magnitude of neutrino masses that is 
suggested by the discovery of atmospheric and solar neutrino oscillations. One of the most crucial 
and generic predictions of grand unification, however, is that the proton must ultimately decay into 
leptonic matter such as a positron and a meson, revealing quark-lepton unity. A class of well- 
motivated theories of grand unification, based on the symmetry of SO(IO) and Supersymmetry, 
which have the virtue that they successfully describe the masses and mixings of all quarks and lep- 
tons including neutrinos, and which also explain the origin of the excess of matter over anti-matter 
through a process called "leptogenesis", provide a conservative (theoretical) upper Umit on the 
proton lifetime which is within a factor of ten of the current lower limit. This makes the discovery 
potential for proton decay in a next-generation experiment high. 

From a broader viewpoint, proton decay, if found, would provide us with a unique window to 
view physics at truly short distances - less than 10~^^ cm., corresponding to energies greater than 
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10 ° GeV - a feature that cannot be achieved by any other means. It would provide the missing 
link of Grand Unification. Last, but not least, it would help ascertain our ideas about the origin 
of an excess of matter over anti-matter that is crucial to the origin of Ufe itself. In this sense, and 
given that the predictions of a well-motivated class of Grand Unified theories for proton lifetime 
are not far above the current limit, the need for an improved search for proton decay through a 
next-generation detector seems compelling. The theoretical guidance provided by some promising 
models points towards the need for improved searches for proton decaying into vK~^ and e'^Tt^ 
modes with lifetimes less than about 2 x 10^^ and 10^^ years, respectively. Should proton decay be 
discovered in these modes, valuable insight would be gained by searches for other related modes 
including fi'^n^ and fi'^K^. 

Current status of experimentation: The "classical" proton decay mode, p — > e'^Tt^, can be ef- 
ficiently detected with low background. At present, the best limit on this mode ( > 5.4 x 10^^ 
yr, 90% CL) comes from a 92 kTon-yr exposure of Super-Kamiokande. The detection efficiency 
of 44% dominated by final-state 7t^ absorption or charge-exchange in the nucleus, and the ex- 
pected background is 2.2 events/Mton-yr. The mode p — > v^+, is experimentally more difficult 
in water Cherenkov detectors due to the unobservable neutrino and the fact that the kaon is below 
Cherenkov threshold. The present limit from Super-Kamiokande is the result of combining several 
channels, the most sensitive of which is K'^ — > jU+ V accompanied by a de-excitation signature from 
the remnant ^^A^ nucleus. Monte Carlo studies suggest that this mode should remain background 
free for the foreseeable future. The present Umit on this mode is > 2.2 x 10^^ yr (90% CL). 

Requirements for the next stage of experimentation: Since the lifetime of the nucleon is un- 
known, and could range from just above present limits to many orders of magnitude greater, 
increases in sensitivity by factors of a few are insufficient to motivate new experiments. Thus, 
continued progress in the search for nucleon decay inevitably requires much larger detectors than 
Super-Kamiokande. The efficiency for detection of the e'^Tt^ mode is dominated by pion absorp- 
tion effects in the nucleus, and cannot be improved significantly. An order of magnitude improve- 
ment in this mode can only be achieved by running Super-Kamiokande for an additional 30-40 
more years, or by constructing an order of magnitude larger experiment. The decay modes of the 
nucleon are also unknown, and produce quite different experimental signatures, so future detec- 
tors must be sensitive to most or all of the kinematically allowed channels. Moreover, the enor- 
mous mass and exposure required to improve significantly on existing limits (and the unknowable 
prospects for positive detection) underline the importance of any future experiment's ability to 
address other important physics questions while waiting for the proton to decay. 

New facilities under consideration: A variety of technologies for discovery of nucleon decay 
have been discussed. Of these, underground water Cherenkov appears to be the only one capable 
of reaching lifetimes of 10^^ years or greater. Cooperative, parallel studies of a future underground 
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water Cherenkov proton decay experiment are underway in the U.S. and Japan. The proposed de- 
signs range from 300 kTon (14 times Super-Kamiokande) to 1 Mton. Liquid Argon or scintillation 
techniques have also been discussed in the proton decay community and may have significant effi- 
ciency advantages for certain modes that are dominant in a certain broad class of SUSY theories. 
Liquid Argon time projection chambers potentially offer very detailed measurements of particle 
physics events with superb resolution and particle identification. Liquid Argon feasibility will be 
demonstrated in the near future with the operation of a 600-ton ICARUS detector. If expectations 
are correct, it should have a sensitivity that is equivalent to a 6000-ton water Cherenkov detec- 
tor in the p — > v^"*" mode. The liquid scintillator approach is presently being explored with the 
1 kTon KamLAND experiment. It should also have enhanced sensitivity to this mode by directly 
observing the by dE/dx and observing the subsequent — > jU"'"v decay. 

Performance and feasibility: Detailed Monte Carlo studies, including full reconstruction of sim- 
ulated data, indicate that the water detectors could reach the goal of an order of magnitude improve- 
ment on anticipated nucleon decay limits from Super-Kamiokande. With sufficient exposure, clear 
discovery of nucleon decay into would be possible even at lifetimes of (few) x 10^^ years 
where present analyses would be background-limited, by tightening the selection criteria. For in- 
stance, with a detection efficiency of 18%, the expected background is only 0.15 events/Mton-yr, 
ensuring a signal:noise of 4:1 even for a proton lifetime of 10^^ years. A water Cherenkov detector 
would also provide a decisive test of super-symmetric SO(IO) grand unified theory by reaching a 
sensitivity of a (few) x 10^"^ years for the vK'^ mode. 

As we have discussed, a much smaller liquid argon could do particularly well on the mode vK^ 
as the efficiency could be as much as 10 times larger than that in the water Cherenkov detectors 
due to the extraordinary bubble chamber-like pattern recognition capabilities. Due to this, a single 
observed event could be powerful evidence for a discovery. The e'^Tt^ mode however would be 
limited by the smaller size of these detectors. 

The search for n-nbar oscillation is another test of baryon non-conservation. While this is not 
one of the favorite predictions of conventional SUSY grand unification, this process, taking place 
in the nuclear potential, can reach an equivalent sensitivity to baryon non-conservation of 10^^ 
years. 

14.2. Observation of Natural Sources of Neutrinos 

All of the detector technologies we consider will lead to enhanced detection and study of neutri- 
nos from natural sources such as the Sun, Earth's atmosphere and lithosphere, and past and current 
supernova explosions. There may also be previously unsuspected, natural neutrino sources that 
appear when the detector mass reaches the hundreds of kilotons scale. The liquid scintillator tech- 
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nique is of particular note here because it could allow the detection of low energy antineutrinos 
from Earth's lithosphere. This physics, however, requires low energy thresholds which are diffi- 
cult to obtain without eliminating cosmogenic background by locating the detector at great depth 
and with careful selection of materials with low intrinsic radioactivity for the detector construction. 
The low activity concerns become important if we attempt to push the threshold to below ~ 5MeV. 
The low activity requirement is not essential for accelerator physics. 

Solar neutrinos have already been observed in the Super- Kamiokande and SNO detectors. If the 
large detector concepts discussed here result in construction of the underground experiment, it may 
become possible to increase the observable event rate enough to clearly observe spectral distortion 
in the < 5 to 14 MeV region. One could also measure the as yet undetected hep solar neutrinos 
(with an endpoint of 18.8 MeV) well beyond the endpoint ( 14 MeV). These measurements 
would require a very comprehensive understanding of the detector systematics and energy resolu- 
tion, but a better determination of the solar spectrum as well as detection of the day-night effect 
with high statistics would represent a significant advance in the evolution of solar nuclear physics 
measurements. 

The observation of supernova neutrino events in a large neutrino detector of the type being 
discussed in this report is straightforward and has historical precedent. The SN 1987 A super- 
nova, in fact, was seen by two large water Cherenkov detectors (11 events in Kamiokande-II (total 
mass 3kT) and 8 events in 1MB (total mass 7kT)) that were active in proton decay searches at that 
time. The predicted occurrence rate for neutrino- observable supernovas (from our own galaxy and 
of order 10 kpc distant) is about 1 per 20 years, so events will be very rare. However, the informa- 
tion from a single event, incorporating measured energies and time sequence for tens of thousands 
of neutrino interactions, obtained by a very large neutrino detector, could provide significantly 
more information than has ever been obtained before about the time evolution of a supernova. 
In addition to obtaining information about supernova processes, the small numbers of SN 1987a 
neutrino events have been extensively used to limit fundamental neutrino properties. Supernova 
processes continue to have very high interest because of the recent detection of the acceleration of 
the rate of expansion of the universe using type la supernova. Recent work has shown that diffuse 
neutrino events from past core collapse supernova (which produce neutrino bursts) could be used 
to gain independent knowledge on the cosmological evolution parameters [|79l . Therefore detection 
of supernova neutrinos, either as a burst from a single supernova or as a diffuse source from past 
supernovas, should be a key mission of the multipurpose detector facility. 

Recently, there has been substantial progress in the detection of relic supernova neutrinos using 
the inverse beta decay reaction Vg + p ^ e+ + n. This ability could be obtained with some futher 
investment as described below. Neutrons are presently invisible in water Cherenkov detectors. Af- 
ter thermalizing, they are captured by free protons in the water, emitting a 2.2 MeV gamma which 
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is well below a typical threshold and which is also overwhelmed by the large radon backgrounds at 
such energies. However, if we were to dissolve gadolinium in the form of gadolinium (tri)chloride, 
GdClTi, in the water (the price of gadolinium has dropped three orders of magnitude in recent years, 
making such a detector affordable) then the experiment would become sensitive to the neutron cap- 
ture gamma cascade (total energy = 8.0 MeV) produced by Gd following positron emission from 
the inverse beta reaction llSOll . With a concentration of 0.1% Gd (0.2% GJC/3) by mass, over 90% 
of the neutrons will be visibly captured on Gd rather than on protons. 

By requiring coincident signals, i.e., a positron's Cherenkov light followed shortly thereafter 
(< 100/i5) and very close to the same spot by the gamma cascade of a captured neutron, back- 
grounds to the diffuse supernova neutrino signal could be greatly reduced. Diffuse supernova neu- 
trino background (DSNB) models vary, but with the Gd in the water the SOkton Super- Kamiokande 
should see about five DSNB events each year above 10 MeV with essentially no background. One 
can easily imagine a next-generation water Cherenkov detector seeing > 100 supernova relic neu- 
trinos every year. Adding gadolinium would greatly improve the response to a supernova within 
our own galaxy as well, allowing the deconvolution of the various neutrino signals (charged cur- 
rent, neutral current, elastic scattering) and, among other things, doubling the pointing accuracy 
back to the progenitor star. Such a detector would also be sensitive to late black hole formation 
to much longer times than at present, since the distinctive coincident inverse beta signals can be 
distinguished from the usual singles backgrounds. An abrupt cutoff of these coincident signals 
would be the unmistakable signature of a singularity being born. 

The continued study of atmospheric neutrinos in the large underground detector will provide 
useful additions to the program carried out so successfully by the Super- Kamiokande Experiment. 
A detector with mass (~1 Mton) would be a powerful tool for studying neutrino physics from 
atmospheric neutrinos. Thanks to the larger dimensions of the detector, higher energy neutrino- 
induced muons can be fully contained and their energy can be measured. Using the atmospheric 
neutrino flux, the distinctive oscillatory pattern as a function of L/E could be directly observed. The 
factor of 20 increase in detector fiducial mass will allow statistical improvements in all the topics 
studied and, perhaps, the emergence of new scientific topics. Other natural sources of neutrinos, 
such as lithospheric neutrinos, have not yet been studied extensively and could, in principle, be 
observed by the new detector concepts. An initial result in this area has recently been announced 
by KamLAND. Typically, the neutrino energies for these processes are below 10 MeV and are 
sensitively dependent upon the low-energy threshold capability of the new detectors. The liquid 
scintillator detector concepts are likely to have the best opportunities for advancing these topics, 
but liquid argon detectors could also contribute. 

Finally, we note that there may be galactic sources of neutrinos that are of lower energy and 
greater abundance than the ultra high-energy neutrino sources to be explored by detectors such 
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as the Ice Cube Cherenkov detector now being constructed deep under the Antarctic ice sheet 
by an NSF sponsored collaboration. Galactic neutrinos have a natural source in inelastic nuclear 
collisions through the leptonic decays of charged secondary pions. This source is expected to be of 
comparable intensity and energy distribution to the high-energy photons that are bom from neutral 
pion decays in the same collisions. Such neutrino sources, currently not detectable with Super- 
Kamiokande, could be seen by a megaton-class neutrino detector that runs for several decades. 

14.3. Depth requirements for non-accelerator physics 

It is difficult to consider all possible non-accelerator physics channels and precisely predict 
the most optimum depth for either water Cherenkov detector or a liquid argon detector. The an- 
swer could easily depend on various technical assumptions, but it is certainly clear that depth 
comparable to or larger than present detectors (Super- Kamiokande is at 1000 m of rock or 2700 
meter-water-equivalent depth) is needed for the best physics reach. A quantitative summary of 
depth considerations can be seen in [fSTIl . 

Nucleon decay modes can be divided in two classes: ones where all of the nucleon energy is 
visible and ones where some of the nucleon energy escapes detection. In the first case, the total 
momentum and energy balance is a powerful tool for background reduction, and it has been often 
argued that these modes should require only modest shielding from cosmic rays. Indeed, most of 
the decay modes that were searched for in the first generation detectors required only modest depth. 
1MB operated successfully at a depth of 2000 feet. However, in a very large water Cherenkov 
detector, cosmics not only produce background, but also reduce the live-time of the experiment 
by keeping the detector occupied by frequent large energy deposits. If we require live-time to be 
more than 90%, a shallow depth of few tens of meters appears sufficient. This conclusion does not 
include consideration of the data rate, which is continuous for non- accelerator physics, and could 
be unmanageably high near the surface. The requirement of a reasonable data rate (< lOHz of 
muons) increases the depth required to approximately the Super- Kamiokande depth. 

For the second class of nucleon decays in a water Cherenkov detector, a low energy tag from 
dexcitation photons may need to be used (For example p vK^ with a ~ 6.3 MeV gamma from 
^^A'^ de-excitation followed by jl^v with lifetime of 12 ns). These require low energy 

thresholds for photons. This is difficult with a background of fast-neutron (spallation products 
from muons in the detector or in the surrounding rock) induced low energy background events at 
shallow depths. Nevertheless, since the tagging photon is in-time to the main event (with time 
window of < 50ns), one could conclude that these events also may not require much more than 
Super-Kamiokande depths. A subclass of events are, however, subject to fast neutron backgrounds. 
As an example of this, the mode n vvv can be searched for by observing the de-excitation of 
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the residual nucleus. The proposed ultimate DUSEL depth (about 6500 mwe) would reduce the 
muon background by about a factor of 100 with respect to Super- Kamiokande and certainly help 
in the observation of these modes with a low energy component. 

For a liquid argon calorimeter, much higher resolution may permit relaxation of these issues. 
In particular, the vK~^ mode could be much easier to detect because the kaon could be identified 
by its energy deposit (dE/dx). Nevertheless, some minimum depth will very likely be necessary to 
reduce backgrounds from fast neutrons and to reduce the data rate to manageable levels. 

For solar neutrinos in a water Cherenkov detector, the important issue is dead-time introduced 
by spallation induced fast neutron backgrounds. At Super-Kamiokande this dead-time is ~ 20%. 
To maintain the same level of dead-time for a much larger detector, depth similar to or greater than 
Super-Kamiokande (2700 mwe) will be needed. For a liquid argon detector, this requirement could 
be relaxed because the dead volume around a cosmic muon could be better defined. 

For a supernova in our galaxy (lOkpc), the signal level is so large (~ 10000/ sec over a 10 
sec burst), that the spallation background at depths as shallow as 500 mwe are manageable. For 
detection of supernova in neighboring Andromeda (~750 kpc), however, greater depth (> 1300 
mwe) is needed. Optimizing depth for diffuse relic supernova neutrino search needs to take into 
account the deadtime loss as well as background from spallation products such as ^Li which beta 
decays and then ejects a neutron. The analysis in [[811 suggests that this search may require depths 
similar to Super-Kamiokande even if one could get the enhancement in signal to background from 
gadolinium loading. 

In summary, the driving issues for depth consideration for future large water Cherenkov or 
liquid argon detectors will be backgrounds to low energy events from spallation products and data 
rates. If one wants to maintain sensitivity to specific important physics channels such as p ^ vK~^ 
in a water detector, and solar and supernova neutrinos in either technology, depth in the same range 
as the current Super-Kamiokande detector is needed. Greater depth will enhance the physics reach 
of the detector. 

15. RESULTS AND CONCLUSIONS 

The following summary results and conclusions were discussed at the Sep. 17 2006 meeting of 
the study group. The broad conclusions have been refined by significant additional numerical work 
since then. We have first listed the broad conclusions from the study. A summary of comparisons 
for the various experimental approaches follows. 

• Very massive detectors with efficient fiducial mass of > 100 kTon (in the case of water 
Cherenkov several hundred kTon and in the case of a liquid argon detector ~ 100 kTon; 
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for accelerator based neutrino physics these two would be roughly equivalent in sensitivity) 
could be key shared research facilities for the future particle, nuclear and astrophysics re- 
search programs. Such a detectors can be used with a long baseline neutrino beam from an 
accelerator laboratory to determine (or bound) leptonic CP violation and measure all param- 
eters of 3 generation neutrino oscillations. At the same time, if located in a low background 
underground environment, they would have additional physics capabilities for proton de- 
cay and continuous observation of natural sources of neutrinos such as supernova or other 
astrophysical sources of neutrinos. 

The Phase-II program will need considerable upgrades to the current accelerator intensity 
from FNAL. Main Injector accelerator intensity upgrade to ~ 700 kW (from the current 
~200 kW) is already planned for Phase-I of the program (NOvA). A further upgrade to 1.2 
MW is under design and discussion as described briefly in this report. The phase-II program 
could be carried out with the these planned upgrades. Any further improvements, perhaps 
with a new intense source of protons, will obviously increase the statistical sensitivity and 
measurement precision. Such an upgrade could significantly reduce the running times (es- 
pecially in antineutrino mode) and increase statistical precision. 

A water Cherenkov detector of multi-lOOkTon size is needed to obtain sufficient statistical 
power to reach good sensitivity to CP violation. This requirement is independent of whether 
one uses the off-axis technique or the broadband technique in which the detector is housed 
in one of the DUSEL sites. 

High signal efficiency at high energies and excellent background reduction in a liquid argon 
TPC allows the size of such a detector to be smaller by a factor of 3 to 4 compared to a water 
Cherenkov detector for equal sensitivity. Such a detector is still quite large. 

The water Cherenkov technology is well-known. The issues of signal extraction and back- 
ground reduction were discussed and documented at length in this study. The needed back- 
ground reduction and energy resolution is achievable and well understood for the broadband 
beam approach, but not yet fully optimized. Key issues for scaling up the current generation 
of water Cherenkov detectors (Super-Kamiokande, SNO, etc.) and locating such detectors 
in underground locations in DUSEL have been investigated. The cost and schedule for such 
a detector could be created with high degree of confidence. A first approximation for this 
was reported to this study. 

For a very large liquid argon time projection detector key technical issues have been iden- 
tified for the building of the detector. A possible development path includes understanding 
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argon purity in large industrial tanks, mechanical and electronics issues associated with long 
wires, and construction of at least one prototype in the mass range of 1 kTon. 



In the course of this study we have examined the surface operation of the proposed massive 
detectors for accelerator neutrino physics. Water Cherenkov detector are suitable for deep 
underground locations only. Surface or near surface operation of liquid argon TPCs is pos- 
sible, but requires that adequate rejection of cosmic rays be demonstrated. Surface or near 
surface operation capability is essential for the off-axis program based on the existing NuMI 
beam-line because of the geographic area through which the beam travels. 

For an off-axis program based on the NuMI beam-line, baselines of about 800 km and off- 
axis distances of 10 to 40 km were considered for CP violation physics. Since the detector 
location is on the surface the best choice appears to be a fine grained detector such as a large 
liquid argon TPC. The scenarios considered for this program were: a) 100 kTon LArTPC at 
the 2nd oscillation maximum (40-60 mrad) in conjunction with the Phase I NOvA detector, 
b) 100 kTon LArTPC at the Phase I NOvA site. For scenario a) we find that the simple 
addition of a 2nd detector does not have significant sensitivity for CP. Scenario b) does have 



sensitivity as shown in Figure 24 A third scenario, using two detectors of 50 kTon each at 



the first and second maximum has also been analysed (see Figure |26]). Additional scenarios 
are presented in Table |X} 

For a wideband program to DUSEL (either at Henderson (1495 km) or Homestake(1290 
km)), two choices for detector technology were considered: a deep sited large water 
Cherenkov detector with fiducial mass of ~ 300 kT or a 100 kT liquid argon TPC (which may 
be located either on the surface or underground). These were found to have good sensitivity 
for CP violation after exposure to the same amount of beam. The better signal to background 
ratio for the liquid Argon TPC allows for better sensitivity which can be compensated by in- 
creased exposure or a larger water Cherenkov detector. The sensitivity for 1300 km location 



and its variation for exposure are shown in Figures 11 to 22 The sensitivity was found to be 
about the same for 1495 km. 

Baselines shorter than 500 km on the NuMI beamline from FNAL have severe technical 
limitations for performing the CP violation science because of the low energy of the oscil- 
lated events, difficulty of separating the ambiguities due to mass hierarchy, and the surface 
location of the massive detectors. 
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15.1. Brief comparison of experimental approaches 



In the course of this year long study we have been able to draw several very clear conclusions. 
Regardless of which options evolve into a future program, the following will be required. 

1. A proton source capable of delivering 1 - 2 MW to the neutrino production target. 

2. Neutrino beam devices (targets and focusing horns) capable of efficient operation at high 
intensity. 

3. Neutrino beam enclosures which provide the required level of environmental and personnel 
radiological protection. 

4. Massive (>>100 kton) detectors which have have high efficiency, resolution and back- 
ground rejection. 

5. For each of the above items, significant investment in R and/or D is required and needs to be 
an important aspect of the current program. 

We have found that the main areas of this study can be discussed relatively simply if we divide 
them into two broad categories : 1) The neutrino beam configuration and 2) The detector technol- 
ogy. Further, we are able to summarize our conclusions in two tables which show the pros and 
cons of the various options. 



In Table XII we compare the pros and cons of using the existing NuMI beam and locating 
detectors at various locations, versus a new wide band neutrino beam, from Fermilab but directed 
to a new laboratory located at one of the potential DUSEL sites, i.e. at a baseline of 1300 to 2600 
km. 



In Table XIII we compare the pros and cons of constructing massive detectors ( 100 - 300 kT 
total fiducial mass) using either water Cherenkov or liquid argon technology. 



15.2. Project timescales 

In the following we briefly comment on the possible timelines for the different components of 
the program we have described in the report. At this stage it is difficult to understand the funding, 
manpower, and other constraints to the program, therefore the study group has decided to comment 
only on technically driven schedules. 

• The FNAL proton upgrade timeline: The SNuMI project which aims to upgrade the Fer- 
milab accelerator complex to deliver higher intensity from the Main Injector, submitted a 
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Pro 


Con 


NuMI On-axis 


Beam exists; 
Tunable spectrum; 


L ~ 735 km 

Sensitivity to mass hierarchy is limited 
Difficult to get flux < 3 GeV 


NuMI Off-axis 
(1st maximum) 


Beam exists ; 
Optimized energy; 
Optimized location for 
1st detector; 

Site will exist from NOvA project; 


L ~ 800 km 

Limited sensitivity to mass hierarchy 


NuMI Off-axis 
(2nd maximum) 


Beam exists; 
Optimized energy; 
Improves mass hierarchy 
sensitivity if 613 is large; 


L ~ 700-800 km; 
Extremely low event rate; 
A new site is needed; 
Energy of events is ~ SOOMeV; 
Spectrum is very narrow 


WBB to DUSEL 


More optimum (longer) baseline; 
Can fit oscillation parameters 
using energy spectrum; 
Underground DUSEL site for detector; 
Detector can be multi-purpose; 


New beam construction project >$100M; 
Multi-year beam construction; 



TABLE XII: Comparison of the existing NuMI beam to a possible new wide band low energy (WBLE) beam 
to DUSEL 



conceptual design report (CDR) in the fall of 2006. The timescale for the project will be- 
come clearer after the review process is completed. A preliminary timeline has been pro- 
vided to this study. The complete upgrade will be carried out in two steps. In the first step, 
the recycler based upgrade (proton plan phase-II) will bring the total beam power to 700 kW 
by early 2011. In the second step, the accumulator upgrade (the complete SNuMI project) 
will bring the total intensity to 1.2 MW. An aggressive plan calls for performing the com- 
plete upgrade up to 1.2 MW by 2012. But this will depend on the outcome of reviews and 
discussions that will take place in the next year. 

Construction of a new beam towards DUSEL: Construction of a possible new beam to- 
wards is not part of the SNuMI project. Only preliminary discussions, cost, and schedule 
estimates exists. The scope of the project is similar to the NuMI project which was de- 



scribed in Section 7.1 Based on the NuMI experience, a rough outline for the project could 
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Pro 


Con 


Water 


Well understood and proven technology; 


Must operate underground; 


Cherenkov 


Technique demonstrated by SuperK (50kT); 


Scale up factor is < 10; 






Cavern stability must be assured 






and could add cost uncertainty; 




New background rejection techniques 


NC background depends on spectrum 




available; 


and comparable to instrinsic background; 




Signal energy resolution ~ 10%; 


Low Vf, signal efficiency (15-20%); 




Underground location 






makes it a multi-purpose detector; 






Cosmic ray rate at 5000ft is ~0. 1 Hz. 






Excellent sensitivity to p ^ 7r"e+ 


Low efficiency to p ^ K^v 


Liquid 


Technology demonstrated by 


Scale up factor of ~300 is needed; 


Argon 


ICARUS (0.3kT); 




TPC 




Needs considerable R&D for costing; 




Promises high efficiency and 


Not yet demonstrated by 




background rejection; 


simulation of a large detector; 




Has potential to operate 


Needs detailed safety design for 




on (or near) surface; 


deep location in a cavern; 




Could be placed on surface 


Needs detailed demonstration 




either at NuMI Offaxis or DUSEL; 


of cosmic ray rejection; 






Surface cosmic rate ~500kHz; 




Better sensitivity to 


Surface operation limits 






physics program; 



TABLE XIIL Comparison of Water Cherenkov to Liquid Argon detector technologies 



be: 1 to 2 years for preparation and geological site investigations, 2.5 to 3 years for civil 
construction, and 1 year for installation of technical systems: a total of 4.5 to 6 years for 
construction of the beam-line. There are a number of issues that are different between NuMI 
and a new beam-line to DUSEL. These are related to the greater downwards angle of the 
DUSEL beam-line and the proximity of the DUSEL beam-line to the FNAL site boundaries. 
These issues and their mitigation will be addressed in a separate note [|6T1 . 

• Construction of a deep large water Cherenkov detector: There are two well recognized 
considerations that define the time scale over which a large water Cherenkov detector could 
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be built: the underground cavern construction and manufacturing of large numbers of photo- 
multiplier tubes. For both the single cavern and the multiple cavern concepts of the detector 
a significant period of exploratory excavations and bore holes will be needed. After this 
period (~ 1 — 2 yrs) approximately 5 to 6 yrs of excavation is needed to reach the needed 
total volume. The PMT manufacturing period depends on the choice of the PMT, which is 
different for the two different concepts for the detector. For 20 inch PMTs, the UNO plan 
calls for manufacturing 56000 tubes in about 8-10 yrs. For the Homestake multiple module 
proposal, the plan calls for manufacturing 10-12 inch PMTs with a rate of about 150000 
tubes in 6-7 yrs. The collaborations are communicating with two large manufacturers of 
hemispherical PMTs. Preliminary conclusions are that each of the two manufactures have 
sufficient capacity currently to produce 10-12 inch diameter PMTs at about 1/2 the rate that 
is needed for these projects. For either choice, smaller or larger diameter, the production 
capacity needs to be enhanced to meet the need, but the investment needed in not considered 
extraordinary. There could be bottlenecks in production of materials (for example, glass) 
that need to be fully understood. 

• Construction of a very large LARTPC: The cost and schedule estimate for a very large 
liquid argon TPC of size (50 to 100 kTon) must be preceded by a series of development steps. 
Although the viability of the technique has been established by the ICARUS group, a factor 
of 10 cost reduction is required to make a very large detector economically possible. The 
development program is outlined in [|65l and contains three projects. One project involves 
techniques for the purification of liquid argon to achieve long electron drift times, low noise 
electronics design, and materials qualification. A second project is the construction of a ~3 
ton module to test design concepts for the very large detector; and the third project is the 
design and construction of a ~ 1 kT detector to be constructed using the techniques proposed 
for the very large detector. 

The first project is in progress at FNAL and Yale. Long (many millisecond) electron drift 
lifetimes have been achieved and the project is expected to be complete by mid 2007. Depen- 
dent on funding the second project could produce results by the end of 2007. The siting and 
mass of the detector to be proposed for the third project are under discussion. Once a choice 
is made, the group would like to start the design immediately. Completion of the design for 
project 3 is expected to take 1 year and requires successful completion of the other projects. 
A preliminary cost for project 3 at this time is ~ $10 M. The above program is essential for 
a LARTPC detector on the surface or underground. The cost implications for siting a very 
large detector at underground locations are being discussed, but they need further work. 
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APPENDIX A: ANSWERS TO QUESTIONS RAISED BY NUSAG 



1. Noting the existence of discrepant sensitivity calculations even for the same detector, it 
would be most useful to have any such calculations performed with consistent assumptions 
and methodologies. 

a) Fixed, common, stated values of the mixing parameters not explicitly under study. 

b) Common, stated and plotted, cross sections vs. En. Common, stated nuclear models. 

c) Stated assumptions about energy resolution, background rejection. 

d) if appropriate, common total p.o.t. If sensible, use a common proton energy and anti-nu 
running fraction. If not, state the optima chosen. 

e) What methods are used to extract the oscillation parameters from the final event sample 
(counting? fitting the spectrum?) 

f) Standardized, stated method for defining sensitivity. 

2. Give sufficient detail in tables and/or plots to allow a reader to understand how the 
numbers for rates or sensitivities are obtained. We would expect that many of the results 
would be easily accessible to a physicist with a calculator. Here are some useful inputs that 
come to mind (meant as a guide only): 

a) Specify the signal channel(s). (We will assume here that it is quasi-elastic.) b) What 
simple cuts (energy, etc.), if any, do you apply? 

c) The number of INTRINSIC Vg events reconstructed as signal, and their reconstructed 
energy spectrum (in reconstructed Enu(QE) or Evis, or Ee, or whatever you'll use.) 

d) What is the purity of the QE selection, that is, for true Vg events, what fraction of those 
selected as QE are actually QE (as a function of E)? 

e) The total number of NC 7t^ events, and spectra vs. true Enu and tu^ momentum. 

f) The number of NC 7t^ events reconstructed as signal, and their reconstructed energy 
spectrum. What is the true Enu spectrum for the NC piO events reconstructed as signal? 

g) The NC 7t^ rejection assumed, as a function of... (;r^ momentum?) 

h) The assumed systematic errors on each of the backgrounds, with any relevant depen- 
dence on energy. How are these estimates arrived at? 

i) The assumed signal efficiency as a function of energy. How are these estimates arrived 
at? 

j) Provide tables and spectra (vs. true and reconstructed Ey) giving the initial population 
of events, before cuts, by process (QE, CCpi+, DIS,...), how these numbers diminish as the 
cuts are applied, and in the final sample at the various oscillation parameter test points. An 
entry at the 3-s sensitivity limit would be informative. Scatter-plots of reconstructed vs. true 
En for individual signal and background channels may be informative. 
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3. Specify the level of simulation that goes into your currently-generated sensitivity esti- 
mates. For example: 

a) How is energy resolution treated? Give a plot of the assumed energy resolution (electron 
energy and neutrino energy) vs. energy. 

b) How is the selection of QE events treated? 

c) How is the rejection of piO's modeled? 

We are grateful to the NUSAG committee to provide questions that could be used to guide the 
study. The report was written with the desire to answer these questions. Some of the details that 
these questions ask for are in the supporting documents which can be obtained from the study 
website: http://nwg.phy.bnl.gov/fnal-bnl/. To keep the length of the report minimum we have 
decided not to repeat the material that can be found in the body of the report. 

4. What near detector location/size/technology/performance/cost is assumed/needed to 
achieve the assumed systematic errors? 



In section 10. 1 we have summarized the thoughts on the near detector issues for the Phase-II(B) 
DUSEL based broadband approach. The requirement on systematic error on the background are 
relatively modest (10%). The harshest requirement might be on the energy scale systematic of 1% 
which is needed to achieve the best precision on the atmospheric parameters of Amgj and sin^ 2023- 
The main technical issues for the near detector are the location for its deployment, the deviation 
from \/r^ behavior of the flux due to the close location of the detector, and the high event rate at 
the near site. The study did not look at these problems in detail. Fortunately, there is now rich 
experience on these issues from the running NuMI-MINOS experiment. Most of this experience 
can be applied directly to the future project. 

For Phase-II(A) approach using the NuMI offaxis beam, the near detector requirements have 
not been studied. The location of such a detector could be in the existing tunnel that connects the 
NuMI near detector site to the beam tunnel. The study has not looked at the event rates or potential 
difficulties due to the deviation from both l/r^ behavior and from having a source with a wider 
angular acceptance at the near detector than the far detector. 

5. K possible, for comparison purposes, use the same methodologies to make parallel 
sensitivity estimates for NoVA (single detector) and T2K. What sensitivity for NoVA do you 



calculate for the same number of p.o.t. assumed in question 1? Please see Section 13.2 



6. All sensitivity calculations for off-axis configurations must include events from neutri- 
nos in the high-energy peak from kaon decay. 



The detector performance criteria are in Section 13.2 



7. What detector technologies are still worth pursuing for a 2nd off-axis detector - Liquid 
scintillator? Water Cerenkov? Liquid Argon? Other? 

Over the past several years, three potential detector technologies have been considered for a 
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next generation experiment: liquid scintillator (similar to NOvA), water Cherenkov and a liquid 
argon TPC. Here, we summarize the conclusions which have been made to date in regard to the 
detector technology that would be best suited to the off-axis beam. 

Studies of a massive liquid scintillator detector using the simulations developed for NOv A have 
shown that the backgrounds (mostly neutral current) would be approximately 1 : 1 with the signal 
at the second maximum and this option was not considered further. 

A water Cherenkov detector of the size proposed for DUSEL could give sufficient rate in the 
NuMI beam, though there might again be a question of background rejection. However it has been 
concluded that this size of detector must be sited deep underground to avoid being swamped by 
cosmic ray muons and there is no existing deep site available along the NuMI beam, and so we do 
not consider this a viable option. 

A liquid argon TPC has the advantages of high efficiency and high background rejection for 
neutral current events, using the high spatial resolution. Thus for the same sensitivity in the same 
beam it can be factors of around 3 smaller than a water Cherenkov detector. For the sensitivity 
studies we have assumed liquid argon detector(s) with a total fiducial mass of 100 kTon. 

8. There were several references to the possibiUty of a detector at ~250 km in the NuMI 
beam. Is this being pursued by the Working Group? What are the general properties of this 
approach? 

Shorter baseline lengths for NuMI off-axis detectors have been considered in the literature 
[|34ll . For example, for a baseline of 250 km, the first and second oscillation maxima are at 0.50 
GeV and 0.17 GeV, respectively. There are two reasons for considering shorter baselines: small 
matter effects and larger numbers of events because of the closer distance. This solution, however 
has several difficulties. The main ones are: i) The low energies needed forces us to consider 
large off-axis angles (> AOmrad) where the flux of neutrinos is rather poor and the contamination 
from high energy neutrinos from kaon decay large. This largely negates the advantages of the 
larger flux because of the closer distance. The event rate can be easily obtained from [|T2ll by 
scaling, ii) Natural choice for a detector at these energies is a water Cherenkov counter. Since 
most of the events at these energies are quasi-elastics for which a water Cherenkov detector has 
good efficiency, little is gained by utilizing a liquid argon TPC. The water Cherenkov detectors 



needed are too large for operation on the surface as explained in Section 12 iii) For the first 
oscillation maximum, an experiment with almost identical parameters is already being carried out 
in Japan (T2K). Combining the results of T2K and Phase-I of the US program is a subject of 
various reviews [TTTII . 

9. Provide cost and schedule estimates for the same fiducial mass and PMT cover- 
age/channel count used for sensitivity estimates. (We reafize that fiducial/total mass ratios 
may be hard to estimate, but the assumptions should be stated.) 
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We hope that the committee understands that the work reported in this study was carried out in 
parallel in a very short period of time. In addition, members of the study group are considering 
several options for detector sites and design. Therefore it is difficult to obtain complete consis- 
tency in the assumptions that went into simulations versus detector design and cost estimates, etc. 
Obviously we will do the best we can to point out the various points of departure and will depend 
on good judgment. 

The design and cost for a detector in the Henderson laboratory were provided in the presentation 
of Prof. Chang kee Jung at [58]. The fiducial volume for UNO was quoted to be 440 kT at a 
preliminary cost of $437M. 

A conceptual detector design for 300 kT for Homestake was presented to the committee in (59] . 
The authors of that report provided the following answer for their choices: 

A single 100 kiloton module will have a cylindrical fiducial volume with a diameter of 50 me- 
ters and a height of 50 meters. The PMTs on the vertical face of the cylinder will have their 
photo-cathodes on the surface of a 52 meter diameter cylinder. The top and bottom PMTs will 
be separated by 52 meters. This layout defines a fiducial volume that begins 1 meter inside the 
PMT photo-cathode surface. In addition, there will be 0.5 meter veto region surrounding the entire 
detector so that the chamber walls will be on a 53 meter diameter cylinder. 

Our budget estimate for the excavation of the detector chamber was based on a 50 meter diam- 
eter by 50 meter high cylinder. The change from 50 meter to 53 meters involves a volume increase 
of 18% and a surface area increase of 12%. Although our budget breakdown details permit us to 
apply the above scale factors to each of the volume and surface area budget items, we decided, for 
this answer to merely use an average cost increase of 15%. When applied to a single module, the 
construction cost increases from $29.1 million to $33.5 million. Note, that these numbers include 
a contingency of 30%. The total single 100 kiloton detector cost increases from $1 16.6 million to 
$121 million, an increase of 3.6%. 

Similarly, when this cost increase factor is applied to three detectors, the three chamber cost 
increases from $66.1 million to $76 million and the total three detector cost increases from $308.9 
million to $318.8 million. 

In the above we not included the effect of moving the PMTs from the original 50 meter diameter 
cylinder to a 52 meter diameter cylinder, a surface area increase of 8%. If apply this factor to the 
previously assumed PMT and associated electronics cost of $62.1 million this creates another $5 
million increase per 100 kiloton detector. The final cost including all contingencies is then $126 
million for a single 100 kiloton detector and $323.8 million for three such detectors. The above 
increase is less than 10% for budget that has a contingency of about 34%. 

Finally, the simulations for the background estimates were reported in [|38l . They were per- 
formed with the exact geometry of the Super-Kamiokande detector (with 40% PMT coverage 
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using 20 inch diameter tubes). The Homestake detector cost is for 11 inch tubes and 25% cover- 
age. We are confident that this coverage is sufficient because of several factors. First, the PMT 
information both Hamamatsu and Photonis shows that smaller diameter semi-hemispheric tubes 
have higher quantum efficiency (QE) than the 20 inch tube (for example, the Hamamatsu 10.5 inch 
tube has QE of 25% and the 20 inch tubes has QE of 20% at 390 nm) This difference is appar- 
ently well-known and documented. Secondly, the collection efficiency (efficiency of collecting 
the photo-electron into the dynode structure) is also known to be larger for the smaller diameter 
tubes. The collection efficiency factor (an increase of about 13% for the 10.5 inch tubes versus the 
20 inch) is not well documented. Therefore, if corrected by these two effects, the 25% coverage 
with smaller 10.5 inch tubes corresponds to 35% coverage (~ 25% x x 1.13) with the 20 inch 
photomultiplier tubes. We also expect that with a larger detector and far larger granularity, the 
background rejection will get better requiring less total coverage. Nevertheless, we understand 
that all of the above has to be demonstrated with benchtop measurements and detailed simulations 
for which we would like to ask for substantial R&D funds. If we must increase the coverage to 
40% to achieve the physics goals then the cost increase will be approximately $112M which is 
certainly beyond the contingency we have allowed at this point. 

10. For the modular water Cherenkov approach, are you defining 3 modules as your 
basefine detector? 

The authors of report [[59ll reply: 

"Yes, there are three main reasons we believe 3 modules is an optimum choice to start with. 
First, because of the long running times possible at FNAL it appears that a 300 kTon fiducial 
mass is sufficient to reach the desired sensitivity for neutrino oscillations. Second, it is clear 
that for proton decay searches a larger detector is needed, but for current background projections 
a few background events are expected in favored decay modes after exposure of 1 MT-yr. We 
believe that proton decay searches will benefit from further detector and analysis improvements 
after reaching this level of sensitivity. Any modules built after the first 3 modules will benefit from 
this knowledge. Third, there is considerable cost saving by starting the simultaneous construction 
of 3 cavities in the region of relatively well-known Homestake rock near the Ray Davis Chlorine 
chamber as explained in [59]." 

11. For the water Cherenkov counters, we will be eager to hear of progress in algorithms 
for rejecting ;r*^s (and the testing of them). What is the increase in Tt'^ rejection over that 
achieved by Super-K (as a function of ;r° energy) assumed in your current calculations? 
What have you reached with your own simulations/algorithms? Describe briefly the algo- 
rithmic improvements. Does this rejection depend more on total photo-cathode coverage, or 
on granularity? 

The detailed account of the rejection is described in the accompanying report by C. Yanag- 
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isawa et al. (39) as well as in the presentations by Yanagisawa and Dufour [|38l l40l . It is also 
summarized in section [TOTl 

The improvement to the signal to background depends on the neutrino spectrum and neutrino 
oscillation parameters. In the following we use the same neutrino energy spectrum and the same 
neutrino oscillation parameters as used in the above report with the CP violating phase of 45 
degrees. 

For the baseline of 1480 km (Fermilab to Henderson), using the new algorithm the signal to 
the background ratio can be improved from 0.30 to 1.9, while retaining 40% of the signal events 
accepted by the current Super- Kamiokande algorithm. For the baseline of 2540 km (BNL to Home- 
stake), using the new algorithm the signal to the background ratio can be improved from 0.35 to 
2. 1, while retaining 40% of the signal events accepted by the current Super- Kamiokande algorithm. 

Dependence of these results above on the granuarity and photocathode coverage has not been 
studied in a systematic fashion, as we extensively used Super-Kamiokande-I (photocathode cov- 
erage of 40%) Monte Carlo sample. It is also found that, given the 40% photocathode coverage, 
the signal to background ratio can be significantly improved for neutrino events with reconstructed 
neutrino energy of < 1 .2 GeV for a detector with better granularity. Other observations concerning 



this issue are touched upon in section 10.1 



12. Though the worldwide community of proponents of large water Cerenkov detectors 
seems to cooperate in simulations, algorithms, etc., we do not see evidence that there is any 
global planning (site-independent design studies or physics programs, etc.) underway for 
such a detector. Please comment. 

One of the most useful results of the NUSAG process has been the cooperation in simulations 
and algorithms for large water Cherenkov detectors. This cooperation was most evident in the 
participation from the T2KK group in our discussions. We have also had fruitful interactions with 
the proponents of the Frejus based water Cherenkov detectors. 

The description and calculations for the water Cherenkov approach in this report was a result 
of cooperation between two US based groups: the UNO group that wants to develop a single very 
large cavern for the detector and the Homestake based group which wants to develop the detector in 
multiple modules. Both groups have worked together to understand and suppress the backgrounds 
in the detector and also have settled on a similar physics strategy for addressing CP violation in 
neutrino oscillations. 

There are currently 5 well considered proposals for a very large water Cherenkov detector 
worldwide: Hyper-kamiokande detector in Japan, a possible detector in Korea on the same neutrino 
beamline as JPARC to Super-Kamiokande, a very large detector in the Frejus laboratory in France, 
and the two possible sites for DUSEL (Homestake or Henderson mines) in the U.S with either 
a large single volume detector(UNO) or a detector in multiple modules such as the Homestake 
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proposal. A long baseline neutrino oscillation program with emphasis on reaching sensitivity to 
CP violation in neutrino mixing is central to all these proposals. Therefore we believe there is good 
cooperation and agreement on the issue of the physics program for such a detector. 

The other two factors for such a detector are a) site development, and b) photo-sensor and elec- 
tronics R&D and acquisition. The site development is a very large part of this detector design, 
and therefore must be handled locally. There is cooperation and communication between these 
groups to compare costs and schedule for the site development. The costs and schedules appear 
understandable after considering the differences between the engineering and accounting practices 
in these geographic region, but we do not see how global site independent planning can be per- 
formed here. The photo-sensor and electronics R&D is the dominant item in these projects. For 
the photo-sensor R&D we agree that good cooperation could be helpful and lower the costs for 
everyone. 

It should be remarked that each of the above geographical regions has a unique virtue for lo- 
cating this massive detector. For Frejus, it is the availability of CERN as a neutrino source and 
the deep location next to the Frejus highway tunnel. For both HyperK and the Korea based de- 
tectors the uniqueness lies in the location on an existing neutrino beamline from JPARC. For the 
US sites the uniqueness is in first the distance available from Fermilab or BNL (> 1000km) which 
is now recognized as essential for performing the next generation experiment with large CP and 
matter effects, and second the depth available at the potential DUSEL sites to suppress cosmo- 
genic backgrounds. Finally, the size of the detector projects are large but at a scale that could be 
contemplated on a national level. Therefore, global planning for a single such detector and site 
independent studies (in the manner of a very large accelerator project), is perhaps not warranted. 

For Liquid Argon: 

These questions were answered by the liquid argon subgroup. The answers were coordinated 
by Prof. Bonnie Fleming. 

NuSAG recommends that the Liquid Argon group reweight its emphasis from sensitiv- 
ity/reconstruction/pattern recognition to hardware issues and cost estimates. We realize that 
a full switch cannot occur if the LAr group is a big part of the more generic off-axis calcu- 
lations in the Working Group, but, for example, LAr-specific reconstruction and particle ID 
algorithms seem less pressing than technical feasibility. 

13. What has actually been measured on purity of the Ar in a tank made with industrial 
technology? ff not yet tried, when will the first tests be? 

Response: No tests have yet been performed on purity of Argon in a tank made with the in- 
dustrial technology necessary for construction of a massive detector. This test will require a large 
tank, ~1 kTon, constructed using the same techniques as a large detector. This project has been 
envisaged by the LArTPC group as outlined in their report to NuSAG in 2005. A specific plan 
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for this component of the R&D path is presently under study and expected to converge within a 
year. In the meantime, small scale tests using the Materials Test Stand at Fermilab (see writeup for 
details) will have first results addressing purity issues within this year. 

14. When do you expect to have tried 3-m drifts and long wires in the US? What effect 
will the capacitance of very long wires have on electronic noise? 

Response: A program to study 5m drifts using a prototype vessel at Fermilab is in the design 
stages. Depending on funding, results from this project are expected within the next two years. As 
well, long drift tests are underway in Europe on the same timescale. 

A 30m long wire with 4 meters of interconnecting cable to electronics will have a capacitance 
of 620pF [|52l . Using commercial amplifiers, a signal to noise of ~9 can be achieved, adequate 
for LArTPCs. Another configuration that has been considered is to use cold electronics, elimi- 
nating the interconnecting cable. This option is under study at Michigan State University in Carl 
Bromberg's group. 

15. What are the R&D milestones, with an estimated schedule, that would lead to a first 
realistic cost estimate for a detector of the 2nd-off-axis or wide-band class? 

Response: Before developing a realistic cost estimate for a massive detector, 50-lOOkTons in 
size, a reasonably sized, scaled down version of the massive detector should be constructed and 
operated. This detector will test purity in a vessel constructed using the same industrial techniques 
envisaged for the large detector, electronics, ability to handle cosmic ray rate, and cellular design. 
As well, smaller scale tests such as the 5m drift test, long wires test etc, as described in the summary 
document, are necessary. However, it is the IkTon scale test that drives the schedule. The schedule 
for this project is not yet fully fleshed out. 
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APPENDIX B: NUSAG CHARGE 

The charge letter is reproduced on the next two pages 
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Please obtain the letter from http://nwg.phy.bnl.gov/fiial-bnl. 
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Please obtain the letter from http://nwg.phy.bnl.gov/fiial-bnl. 
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APPENDIX C: CHARGE TO TfflS WORKING GROUP 



April 5, 2006 
Dear Colleague, 

This letter is being sent to you as a follow-up to the Long Baseline 
Workshop held at Fermilab on March 6-7. This mailing list is composed 
of those who attended the study and signed up to receive further 
information or have subsequently expressed interest in the 
study. Since the kick off meeting we have redrafted the goals of the 
study. We have inserted a time scale which we judge to be 
achievable. The is goal is described in the attached document. You can 
anticipate that within days you will get a further document in which 
Milind Diwan and Regina Rameika have attempted to parse the study goals 
into a set of work packages . We would like to hear from people who are 
prepared to do some work on these issues . Especially we would be very 
happy to hear from people new to these studies. 

However, as you might expect we do have some likely suspects in mind 
and Gina and Milind will be contacting people to help. Finally, we 
will also be recruiting an Organising/ Advisory Committee to help us 
guide this study. We look forward to seeing progress on this study and 
would welcome your suggestions for additions, adjustments and 
approach . 

With Best Regards, 
Sally & Mont 

Sally Dawson, Chair, Physics Department, Brookhaven National Laboratory. 
Hugh Montgomery, Associate Director, Fermi Natinal Accelerator Laboratory 



100 



U.S. Long Baseline Neutrino Experiment Study 

AprU3, 2006 

Preamble: 

The Fermilab Proton Driver (FPD) Study conducted in 2004/05 concluded that Neutrino Oscillation physics provides 
the main motivation for the FPD, and that the physics case is strong. In addition, the Fermilab 8 GeV linac beam could 
support other physics experiments in parallel with an upgraded NuMI-based neutrino program. The study results were 
presented to the Fermilab PAC in the 2005 Aspen meeting, and were well received. However, the PAC expressed 
interest in the possibility of a fiarther generation of neutrino oscillation experiments at a FPD (beyond NOvA) but 
anchored by the NuMI facility. 

Brookhaven National Laboratory has considered in some detail a very long baseline neutrino oscillations (VLBNO) 
concept using an on-axis, wide band beam but with a very large detector at a longer distance. This beam could 
originate from cither Fermilab or BNL. This approach requires a large underground detector presumably located at the 
NSF's plaimed DUSEL facility, which would also have potential for other fi^ontier physics in addition to neutrino 
oscillation physics. 

While these two approaches have a common goal of understanding neutrino masses and mixings, they arc clearly 
different. Wc would like to have a thorough study and exploration of the differences and potential of the two 
approaches. To that end wc have drafted a charge for a joint Fermilab/BNL snady, the results of which could fom the 
basis for a national program in neutrino physics. This charge follows an initial short workshop which laid out some of 
the issues: 

http://www.ftial.gov/drrectorate/DirReviews/Neutrino_Wrkshp.html 
Timescale: 

The United States neutrino community is heavily engaged in operation and analysis of its existing program. On the 
other hand there are active discussions within advisory bodies and the agencies with a view to setting directions for 
future facilities inside the next year. 

It would be desirable to see results of this U.S. Long Baseline Neutrino Experiment Study before October 2006, 
with a preliminary report by July 15, 2006. 



U.S. Long Baseline Neutrino Experiment Study 

Compare the neutrino oscillation physics potential of 

1 . A broad-band proposal using a either an upgraded beam of around 1 MW from the current Fermilab accelerator 
complex or a fliture Fermilab Proton Driver neutrino beam aimed at a DUSEL-based detector. Compare these 
results with those previously obtained for a high intensity beam from BNL to DUSEL. 

2. Off-Axis next generation options using a 1-2 MW neutrino beam from Fermilab and a liquid argon detector at 
either DUSEL or as a second detector for the Nova experiment. 

Considerations of each should include: 

i) As a fimction of 813, the ability to establish a finite 8 13, determine the mass hierarchy, and search for CP violation 

and, for each measurement, the limiting systematic uncertainties. 

ii) The precision with which each of the oscillation parameters can be measured and the ability to therefore discriminate 

between neutrino mass models. 

iii) Experiment Design Concepts including: 

Optimum proton beam energy 
Optimum geometries 
Detector Technology 
Cost Guesstimate 
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APPENDIX D: STUDY GROUP MEMBERSHIP 



Chairs 

Sarah Dawson (co-chair) 

Hugh Montgomery (co-chair) 



Brookhaven National Lab. dawson@bnl.gov 
Fermi National Accelerator Lab. mont@fnal.gov 



International Advisory Group 



Milind Diwan (co-leader) 


Brookhaven National Lab. 


diwan@bnl.gov 


Regina Rameika (co-leader) 


Fermi National Accelerator Lab. 


rameika@fnal.gov 


Joshua Klein 


University of Texas 


jrk@mail.hep.utexas.edu 


Franco Cervelli 


INFN, Pisa 


franco.cervelli@pi.infn.it 


Maury Goodman 


Argon National Lab. 


maury . goodman @ ANL. GOV 


Bonnie Fleming 


Yale University 


bonnie.fleming@yale.edu 


Karsten Heeger 


Lawrence Berkeley Lab. 


KMHeeger@LBL.GOV 


Steven Parke 


Fermi National Acc. Lab. 


parke@FNAL.GOV 


Takaaki Kajita 


University of Tokyo 


kajita@suketto.icrr.u-tokyo.ac.jp 



A full list of participants is available at http://nwg.phy.bnl.gov/diwan/nwg/fnal-bnl/folks.txtl 
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APPENDIX E: RELEVANT RESOURCES AND URLS FOR THE STUDY GROUP 



Main websites for this study group are: 

http://home.fnal.gov/ rameika/LBL_Study/LBL_mainframe.htm 

http://nwg.phy.bnl.gov/fnal-bnl/ 

Additional materials can be found at: 

http://www.fnal.gov/directorate/DirReviews/Neutrino Wrkshp.html 

http://www.hep.net/nusag_pub/May2006talks.html 

http://www-numi.fnal.gov/ 

http://nwg.phy.bnl.gov/ 

http://www-lartpc.fnal.gov/LBStudy_LAr/2006LB.html 

http://www.dusel.org/ 

http : //www. lbl.gov/nsd/homestake/ 

http://nngroup.physics.sunysb.edu/husep/ 
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APPENDIX F: SCHEDULE OF MEETINGS AND REPORT PREPARATION 



Most of the work of the working group was carried out by small subgroups that worked on the 
individual documents. The work was mostly carried out by email and telephone. The following 
meetings were very helpful for wider interactions. 



November 14, 2005 
March 3, 2006 

March 6-7, 2006 
April 5, 2006 
April 11,2006 
May 20, 2006 
June 27-28, 2006 
July 6, 2006 

July 15, 2006 
September 16-17, 2006 
October 16, 2006 
December, 2006 



FNAL/BNL meeting to explore collaboration, BNL 

Charge letter for NuSAG to examine APS study recommendation for 

a next generation neutrino beam and detector configuration 

First kick-off workshop for organization of the study at FNAL 

Charge letter to the study from Dawson and Montgomery 

Preparation of the task list and assignments 

Presentations to NuSAG committee about the study in Chicago 

Second workshop on detector technologies at FNAL 

Status report to HEPAP from NuSAG, presentation by P. Meyers, 

HEPAP meeting in Washington D.C. 

Deadline for preparation of individual reports from the task list 
Third workshop on preparation of the joint summary report 
Deadline for presentation of the joint report 
Deadline of report from NuSAG to HEPAP 



Other meetings of note where interactions took place are 

June 13-19, 2006 Neutrino 2006, Conference in Santa Fe 

July 11-21, 2006 Neutrino Physics with Liquid Argon TPCs, Yale Univ. 

August 24-30, 2006 NuFact 06 Workshop, UC/Irvine 

Sep 21-23, 2006 NNN06, University of Washington 

March 29-30, 2007 Fermilab Physics Advisory Meeting 
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